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Abstract 

The possibility of structured mixed phases at first order phase tran- 
sitions in neutron stars is re-examined by taking into account charge 
screening and surface effects. The transition from the hadron npe 
phase to the quark phase is studied. Two possibilities, the mixed 
phase and two separate phases given by the double-tangent (Maxwell) 
construction are considered. Inhomogeneous profiles of the electric 
potential and their contribution to the energy are analytically calcu- 
lated. The electric field configurations determine the droplet size and 
the geometry of structures. 



1 Introduction 

It is now commonly accepted that different phase transitions may occur in 
neutron star interiors. The possibilities of pion and kaon condensate states, 
and quark matter state were studied by many authors during the last thirty 
years, see [1, 2, 3, 4] and refs therein. It has been argued that these tran- 
sitions are of first order. They were described in terms of two spatially 
separated phases using the Maxwell construction. Glendenning raised the 
question whether mixed phases in systems composed of charged particles ex- 
ist instead of the configuration described by the Maxwell construction [5] . In 
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particular, the possibilities of hadron - kaon condensate (npe - npeK con d) and 
hadron - quark (H-Q) mixed phases were discussed. The existence of such 
kind of mixed phases in dense neutron star interiors would have important 
consequences for the equation of state, also affecting neutrino emissivities [6], 
glitch phenomena and r modes, cf. [7, 8]. 

Basing on the validity of the Gibbs conditions, in particular on the equal- 
ity of the electron chemical potentials of the phases, refs [5, 9, 10, 8] further 
argued that the Maxwell construction is always unstable. It is due to the 
inequality of the electron chemical potentials of two phases and, thus, due to 
a possibility for particles to fall down from the higher energetic levels char- 
acterized by the higher electron chemical potential of the one phase to the 
lower levels of the other phase. Thereby, refs [5, 9, 10, 11, 8] argued that, if 
first order phase transitions, such as quark, kaon condensate and pion con- 
densate transitions, indeed, occur in neutron star interiors, the existence of 
a wide region of the mixed phase is inevitable. On explicit examples of npe 
- npeK con d and H-Q phase transitions authors demonstrated the energetic 
preference of the mixed phase. 

However, strictly speaking, such an argumentation is in contradiction 
to the results of some other works. It has been recently observed that in 
some models the Gibbs condition of equality of electron chemical potentials 
of two phases can't be fulfilled at all [12, 13, 14], whereas conditions for 
the Maxwell construction are fulfilled. The latter construction is stable in 
these specific cases. A critics of the bulk calculations which ignore finite size 
effects was given in ref. [15] . To include these effects authors used a relaxation 
procedure in which they start with an initial guess for the shape of the electric 
potential in the Wigner-Seitz cell, solve for the kaon equation of motion to 
obtain the charged particle profiles, and then recalculate the electric potential 
using the new profiles. This is then repeated until convergence. Authors 
found that with inclusion of inhomogeneity effects the region of the kaon 
condensate mixed phase is significantly narrowed compared to that obtained 
using standard Gibbs conditions, disregarding finite size effects. Please notice 
that in order the scheme to be completely self-consistent the electric field, as 
the new degree of freedom, should obey the own equation of motion, and it 
enters equations of motions of other fields (including protons and electrons) 
via the gauge shift of the chemical potentials of all charged particles. All the 
charged particle densities are affected by the inhomogeneous electric field 
even in the case, if one artificially suppressed the contribution of the electric 
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field to the charged kaon density. 

For small size droplets, the problem of the construction of the mixed 
phase is analogous to that studied somewhat earlier for matter at sub-nuclear 
densities, cf. [16]. The possibility of the structure, as well as its geometry, 
are determined by competition between the Coulomb energy and the surface 
energy of droplets. Ref. [17] applied these ideas to the description of the 
mixed phase for the H-Q phase transition. The Coulomb plus surface energy 
per droplet of the new phase always has a minimum as a function of the 
droplet radius. On the one hand, this radius should be not too small in order 
for the droplet to have rather large baryon number (A ^> l) 1 and, on the 
other hand, it should be not too large (less than the Debye screening length) 
in order one could ignore screening effects. The Debye screening length for 
the quark matter was evaluated as Ad — 5 fm. The value of the droplet radius 
R m in depends on a poorly known surface tension parameter a. With a small 
value of the surface tension a ~ 10 MeV-fm -2 the droplet size was estimated 
as -Rmm > 3.1 fm and with a ~ 100 MeV-fm -2 , as -R mm > 6.6 fm > Ad- 
Thereby, it was argued that the mixed phase is not permitted in the latter 
case. Corrections to the Coulomb solutions 2 due to screening effects were 
not considered, although one could intuitively expect that for such a narrow 
interval of available values of droplet radii, being of the order of the Debye 
screening length, the screening may significantly affect the results. Also the 
effect of the inhomogeneity of the field profiles in the strong interaction part 
of the energy was disregarded. Therefore, further detailed study of the effects 
of inhomogeneous field profiles, such as screening and surface effects, seems 
to be of prime importance. 

The problem of the construction of the charged density profiles by taking 
into account screening effects is in some sense analogous to that considered 
previously in refs [18, 19, 20] for abnormal pion condensate nuclei. Results 
can be used also for the description of the charge distributions in strangelets 
and kaon condensate droplets. 

The value of the surface tension is poorly known. It has the meaning only 
if there exists a related shortest scale in the problem, being much smaller than 
the typical scale of the change of the electric field. In the case of the kaon 

1 This condition, A » 1, is assumed to be fulfilled. Otherwise quantum effects, like 
shell effects, may affect the consideration. 

2 By the Coulomb solutions we call solutions with the charged density profiles in the 
form of combination of step-functions. Screening effects are disregarded in this case. 
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condensate phase transition there is no necessity to introduce the surface 
tension. One can explicitly solve the equations of motion for all the mean 
meson fields, cf. refs [11, 15]. However one should also include the equation 
of motion for the electric field. The existence of the additional typical scale 
Ad appearing in this equation may affect the solution. In the case of the H-Q 
phase transition there is a natural minimal scale in the problem related to 
the confinement < 0.5 fm. Since the latter scale is much shorter than Ad the 
introduction of the surface tension is well established in the given case. 

With this paper we show that, firstly, the Maxwell construction does not 
contradict the Gibbs condition of equality of the electron chemical poten- 
tials, if one properly incorporates the electric potential. Secondly, inclusion 
of finite size effects, such as screening and surface tension, is crucial for under- 
standing the mixed phase existence and its description. We will consistently 
incorporate these effects. Part of the results of this paper was announced in 
the Letter [21]. 

Our consideration is rather general and main effects do not depend on 
what the concrete system is studied. However we need to consider a concrete 
example to demonstrate the quantitative effects. For that aim we concentrate 
on the description of the H-Q phase transition. Thus we avoid the discussion 
of other possibilities, such as the charged pion condensation, charged kaon 
condensation, see [2, 8] and refs therein. We will consider the simplest case 
of the H-Q phase transition between normal phases. Thus we also avoid the 
consideration of the diquark condensates and diquark condensates together 
with the pion and kaon condensates, see [22] and refs therein. With the 
inclusion of the diquark superfluid gap, the corrections to the effective energy 
functional are expected to be as small as (A/ lib) 2 , where A is the pairing gap 
and ix b is the baryon chemical potential, cf. [14, 23]. However the response 
to the electric field can be a more specific. E.g., in the case of the phase 
transition to the color-flavor-locked (CFL) phase the quark sub-system is 
enforced to be electrically neutral without electrons, cf. [14]. There can be 
even more involved effects, which we do not consider, like the modification 
of the quark and gluon condensates with the increase of the baryon density, 
a modification of masses and widths of all particles, that needs a special 
quantum treatment, etc. 

In sect. 2 we critically review the application of the Gibbs conditions 
to the finite size structures and we discuss the stability of the Maxwell con- 
struction. Since the problem of interpreting the Gibbs conditions and the 
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Maxwell construction is related to the counting of the electric potential from 
different levels we discuss different choices of the gauge in sect. 2. Start- 
ing from sect. 3 we present our results in an arbitrary gauge. In sect. 3 
we develop a general formalism based on the thermodynamic potential. It 
serves as a generating functional, variation of which in the corresponding 
independent variables determines the equations of motion. We demonstrate 
how one can explicitly treat electric field effects. Chemical equilibrium con- 
ditions and the charged particle densities are presented in sect. 4. In sect. 5 
we analytically solve the equation of motion for the electric potential of the 
new phase droplets placed in the Wigner-Seitz lattice and assuming spherical 
geometry we recover corresponding contributions to the energy and thermo- 
dynamic potential. Slab structures are studied in sect. 6. Sect. 7 discusses 
which mixed phase structures, spherical droplets or slabs, are realized in de- 
pendence on the value of the quark fraction volume. The discussion on the 
specifics of the description of the H-Q system is deferred to the Appendices. 
In Appendix A we introduce the model to calculate the volume part of the 
energy of quark and hadron phases. In Appendix B we discuss surface ef- 
fects. In Appendix C we develop a formalism based on the Gibbs potential, 
as generating functional. In Appendix D we study the role of the nonlinear 
correction terms. In Appendix E we discuss peculiarities of the Coulomb 
limit. 

All over the paper we use units H — c — 1. 

2 Gibbs conditions, Maxwell construction 
and electric field effects 

Here, we will illustrate some inconsistencies of a naive treatment of the Gibbs 
conditions. These conditions formulated for spatially homogeneous systems 
were further applied for the description of the Maxwell construction and 
spatially inhomogeneous structures of mixed phase leading to inconsistencies. 

Let us assume that we deal with a first order phase transition in a neu- 
tron star interior between two bulk phases (I) and (II) separated by a thin 
boundary layer. In discussing rather short distance effects one may disre- 
gard spatial changes of the fields due to gravity. In order to describe bulk 
configurations of phases I and II one usually applies the Gibbs conditions 
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[1, 2, 24] 



P l = P 



.a 



(1) 



Pb ~~ 



(2) 



supplemented by the local charge-neutrality relations 



(3) 



Here P a is the pressure, p% is the baryon chemical potential, is the net 
charged density of given phase, a — I, II. 

The configuration corresponding to the coexistence of two locally charge- 
neutral bulk phases is determined by the double tangent construction (Maxwell 
construction). The straight line e(p^) = ^~\i(Pb — Pb) + 6 (Pb) connecting 
the points I and II on the curve of the density dependence of the energy 
density e(p B ) (for temperature T = 0) is characterized by equal derivatives 



and, thus, the pressure equals in both phases. For phase transitions between 
two hadron phases these points relate to the one and the same curve e(p^) 
and correspond to the equal area Maxwell construction on the P{pb) plot. 

Bulk configurations I and II are separated by a boundary layer, which 
presence is ignored in the formulation of the above conditions; e.g., the 
charge-neutrality condition (3) is obviously violated in this layer. Thus, 
the boundary layer is assumed to be rather thin as compared to the sizes 
of regions occupied by both phases. There are two typical sizes characteriz- 
ing the boundary layer between phases. In case of hadron-quark transition 
the shortest scale (ofe) is related to the change of nuclear and quark fields 
in a narrow part of the boundary layer (we call it the surface layer) with 
the length ~ ds < 1 fm and the longest scale (Ad 3> ds) is related to the 
electric field generated in the boundary electric charged layer of the length 
~ A D ~ 5 — 10 fm. Only in the case R 3> A D (R be the minimum of R l 
and -R 11 ), the spatial change of the energy density in this boundary layer 
can be treated as a surface contribution (neglecting corrections of the higher 
order in Xd/R), being rather small compared to the volume one, since the 



d(e/p B ) 
dps 1 



d(e/p B ) 
dps 1 
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surface energy is \d/R times less than the volume energy. Thus, in case 
R ^> Ad, with an accuracy 0(\v/R), one may neglect contributions to ther- 
modynamic quantities of the charged boundary layer and the surface layer 
between two spatially separated bulk phases, as it is usually done in the 
standard formulation of the Maxwell construction. 

The Maxwell construction treatment of first order phase transitions has no 
alternatives for the description of systems with one conserved charge (baryon 
charge in our case). However, neutron stars are composed of charged particles 
and the electric charge is also conserved. Imposing the relations (3) one 
assumes that the electric charge is conserved locally, contrary to the other 
possibility that the electric charge is conserved only globally, see [5, 9]. Thus, 
one should still check a new possibility of the formation of a mixed phase 
constructed of inhomogeneous rather small size charged structures embedded 
in the charge-neutral Wigner-Seitz cells. 

So, the local charge-neutrality conditions might be relaxed, being replaced 
by the global charge-neutrality condition 

Relation (4) introduces the fraction volumes f l of the domain D 1 occupied 
by the phase I and f u , f l = 1 — / n , of the domain D u occupied by the 
phase II. In the case of spherical geometry one deals with spherical droplets 
of the one phase (be phase I) embedded into the other phase (be phase II). 
This is so called the mixed phase. The condition of the total baryon number 
conservation, 

Pb = v- 1 [ dfp B (f) = p\f + p l lf\ (5) 

is also imposed, v = v 1 + v 11 is the total volume occupied by the phases. 

Refs [5, 9] suggested that the conservation of the global electric charge 
requires an extra Gibbs condition for the electron chemical potentials to be 
fulfilled 

A = fi, (6) 

which is formulated in analogy with relation (2). At the first glance, con- 
dition (6) contradicts relations (3), which define p\ ^ /i", as stated in text 
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books. Therefore, refs [9, 10, 8] argued that the double-tangent (Maxwell) 
construction is always unstable since, due to the inequality of electron chem- 
ical potentials, the particles may fall down from the higher energetic level 
characterized by the larger electron chemical potential (for concreteness in 
phase II) to the lower one (in phase I). This argumentation implies the use of 
the same definitions of chemical potentials in Gibbs and Maxwell treatments. 

Now we will show that the Gibbs condition (6) and the Maxwell con- 
dition (3) do not contradict each other since they use different definitions 
of the electron chemical potentials, the global constant quantities fx\ Gibbs 
and /x^Gibbs Gibbs = A^Gibbs = A*e,Gibbs) according to the Gibbs condition 
(6) and the local quantities fj} e loc and /^loc (a4 ioc 7^ A^ioc) according to the 
Maxwell construction, eq. (3). Briefly speaking, the Gibbs condition (6), as 
it was formulated for spatially homogeneous configurations, has no meaning 
in the application to charged systems of a small size, if one does not incor- 
porate electric field effects. It only fixes the level from which one counts the 
electric potential. Therefore, the above mentioned argumentation against 
the Maxwell construction is invalid. It does not take into account effects of 
the electric field arising in any charged systems, at least, near the boundary. 

Thus, first, we need to fix the definition of the electric chemical potential 
of the spatially inhomogeneous system, since we need to recover this quantity 
for phase I based on the knowledge of the quantities of phase II. The two 
ways how one can do it are illustrated in Fig. I. 

The first way (we will call it way I) is as follows. Consider two bulk 
matters occupied by phases I and II separated by a boundary layer (R — \ l D < 
r < R + A"), see Fig.l. The electron chemical potential of the phase I is 
determined as /x e ,Gibbs = /4 Gibbs = A^Gibbs = const in agreement with the 
Gibbs condition (6). The value /^Gibbs * s fixed by the corresponding local 
charge- neutrality condition, see second eq. (3), i.e. /4*Gibbs = /4*ioc = const. 
Then, such a Gibbs condition contains no further information except an 
indication that all the electron energy levels in both phases are occupied up to 
one and the same the top (if we assume that /^g I loc > fi l e loc , where loc is fixed 
by the first local charge- neutrality condition, first eq. (3)) energy level e e = 
/4,Gibbs = /4>c- 111 or der to avoid any contradiction with the local charge- 
neutrality condition for the phase I (first eq. (3)), which would determine the 
value /4 loc 7^ /4,Gibbs> one ne eds to introduce an external electric potential 
well V. The gradient of the latter is necessarily produced in the charged 
boundary layer separating the bulk phases. The value which enters this 
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charge-neutrality condition, second eq. (3), is then V = fi\ loc — ^\oc < 0- ^ n 
the opposite case with //" loc < /ie 5 i oc , the Gibbs condition means that levels 
are counted from the bottom level. The electric potential well is completely 
filled by electrons from the bottom to the top according to the Pauli principle, 
see [18, 19, 20]. Therefore, no instability of such kind arises within the 
Maxwell construction. 

In discussion of the possibility and construction of the finite size struc- 
tures of the mixed phase one also needs to incorporate the electric potential, 
because the small charged droplets have an inhomogeneous profile of the elec- 
tric field. The potential well satisfies the equation of motion, which is the 
Poisson equation, 

AV = 47re 2 ]Tp* h , (7) 

i 

where V = — e0 e is the electric potential well, e is electric charge, e 2 = 1/137, 
and p* h is the charge density of the i particle species. V is the gauge variant 
quantity allowing the replacement V — > V — V° and, thus, being determined 
up to an arbitrary constant V°, see sect. 3. Accordingly, the chemical 
potential is also a gauge dependent quantity. Here in the way I we fixed the 
gauge taking V° — 0, V = V(V = 0), see Fig. 1. Then it is related to the 
electron density p e (r) by means of V(r) = p ei Gibbs _ (37r 2 p e (r)) 1 / 3 . For _R W 
R 3> Ad , where i?w is the radius of the Wigner-Seitz cell, with the Debye 
screening lengths A^ n determined by the same eq. (7), there is no difference 
in solution of (7) for the configuration given by the Maxwell construction and 
for spherical droplets of the mixed phase. What configuration is realized is 
determined by the minimization of the appropriate thermodynamic potential 
over the droplet size. 

The second way (which we will call the way II) is a quasiclassical treat- 
ment also illustrated in Fig. 1. The chemical potential of electrons is intro- 
duced as a local quantity, being unambiguously determined by the value of 
the electric potential well, namely /J> e ,ioc(f) = —V(r). Here we imposed the 
boundary condition V{r) — > —fJ, e , Gibbs at large distances from the droplet, if 
the concentration of droplets f 1 is very small. Thus, the gauge constant is 
taken to be V° = -// e ,Gibbs, i.e. V(r) = V(V° = -// e ,Gibbs) = ~ He,ioc{r) ■ Here 
/-i e ,ioc(j) = (37r 2 p e (r)) 1 / 3 is the quasiclassical quantity, being expressed via the 
local concentration of the particles, cf. [25], p. 321. It is precisely how the 
chemical potential is defined in the usual treatment of Fermi systems. Only 
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in order to distinguish between the local quantity and the global constant 
used in the way I we introduced the subscript "loc" and "Gibbs". Then, 
the local charge-neutrality conditions automatically yield different values of 
/4i oc and //"loc m both phases, as it follows from eq. (3), i.e. in complete 
agreement with the description given by the Maxwell construction. 

The way II chosen, the Gibbs condition (6) should be omitted. One does 
not need to introduce any additional constants, since the value of the top en- 
ergy level is already fixed by the condition V(r) — > — /i e , Gibbs, as r approaches 
the boundary of the Wigner-Seitz cell for the case of a tiny concentration / : 
(i?w ~~ *■ °°)- All the necessary information comes from the Poisson equation 
for the electric potential well, which is the equation of motion for the elec- 
tric field properly derived from the Lagrangian. Equivalently, this equation 
of motion can be obtained by minimizing the appropriate thermodynamic 
potential, being expressed in the corresponding variables. If A^, 11 <C R, there 
exist constant solutions of eq. (7), being valid for each phase outside the 
boundary layer of the length + A". Namely, these constant solutions 
V = const unambiguously guarantee the charge-neutrality conditions (3) 
and, thus, determine the values of the electron chemical potentials of each 
phase, i.e. —V 1 = f/ e loc and —V u = /4*ioc ^ fj} e i oc - I* 1 the boundary layer the 
potential —V(f) or, in other words, the electron chemical potential /x e ,ioc(^0 
varies from the value fj}^ loc to the value f/ e loc . 

An argument that the difference in the values of the electron chemical 
potentials of the two phases produces an instability, since particles may fall 
down from the higher energy level corresponding to the larger value of // ej i oc = 
A^ioc to the lower level (/i e ,ioc — Z-4 ioc) > again does not work. The electric 
field configuration given by the solution of the equation of motion (7) is 
completely filled by the charged fermions from the bottom (V = — a4*i oc ) to 
the top (V = 0) in accordance with the Pauli principle and there is no one 
free state anymore, cf. [18]. If a free state would be formed it would be 
immediately filled by chemical reactions turned out from the equilibrium in 
this case. 

In the case R < A^ the first local charge-neutrality condition (3) is irrele- 
vant. The configuration is determined by the corresponding inhomogeneous 
solution of the Poisson equation. Thus, in order to incorporate the possibility 
of mixed phase structures in the framework of a consistent scheme one needs 
to explicitly introduce the electric field developing according to its equation 
of motion. 
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The equality of the electron chemical potentials of two phases written 
as V 1 = — /4,ioc = ^ H = ~f jl e,ioc would mean that the value of the electric 
potential well V is constant within the volume occupied by the phases. But 
with constant V charged objects can't exist! A charged droplet can be formed 
only, if there is a competition between the electric field energy and the surface 
energy, that permits the existence of a minimum in the Gibbs potential (or 
the thermodynamic potential in the corresponding variables) per droplet at 
some finite droplet size. This should be treated as a necessary condition 
for the existence of a mixed phase. On the other hand, according to eq (7) 
the constancy of the electric potential within some extended space region 
would mean that the r.h.s. of this equation is zero. Therefore it coincides 
with the local charge- neutrality condition, see (3), and we come back to 
the construction with V 1 ^ V 11 describing the case as in the 

Maxwell construction, rather than to a construction of droplets of a small 
size. One definitely should exclude condition (6) in the way II in favor of the 
solution of the equation of motion (7), which then unambiguously determines 
the charged configurations. 

Concluding, the Gibbs condition (6) is automatically satisfied in the way 
I, where the electron chemical potential is constant. This constant then 
enters the r.h.s. of (7). In the way II the electron chemical potential intro- 
duced as Heiocif) is no longer constant and the Gibbs condition (6) should 
be dropped. Both ways are simply related to each other by the gauge trans- 
formation and lead to the very same equations of motion, thus describing the 
same physics (see sect. 3). 

Thus, in the description of mixed phase configurations inhomogeneous 
electric field profiles should be explicitly found. Also surface effects modify 
the Gibbs condition (1). We will discuss the latter effect in Appendix B. 

We consciously paid so much attention to these questions in order to mo- 
tivate the necessity of a consistent account of the electric field effects and to 
remove the so called "contradiction" between descriptions of the mixed phase 
and configurations determined by the Maxwell construction being discussed 
during many years in the literature. 
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3 General formalism. 



Consider the structured mixed phase consisting of two phases I and II. We 
assume droplets of phase I to be located in a lattice described by Wigner- 
Seitz cells. The exterior of the droplets is phase II. Each droplet in the cell 
occupies the domain D l of volume v 1 separated by a narrow surface layer 
D s from matter in phase II ( domain D u of volume v 11 ). We expect particle 
species included in phases I and II to coexist in D s . The thermodynamic 
potential (effective energy) per cell is represented by a density functional 
[26], 

Sl = E\p]-iM] I dfp\-pf [ dffi-iM* I drpl (8) 

JD l JD U Jd s 

E[p] is the energy of the cell 3 , p = {p\, pf} are densities of different particle 
species, i = 1, N 1 in phase I, i — 1, N 11 in phase II, and i — 1, (N l + 
N 11 ) in D$, N 1 , N u are total number of particle species per cell in phases I 
and II. Summation over the repeated Latin indices is implied. We assumed 
that each phase is in the ground state and the matter in phase I or II is 
in chemical equilibrium by means of the weak and strong interactions. The 
term labeled by "S" is the corresponding contribution of a narrow surface 
layer (Ds) separating phases. Equations of motion, Jp- = 0, render 



= ^p « = {I,II,S}. (9) 

The energy of the cell consists of four contributions: 

E[p] = dfe{ in+s M + J m dfe^ n+stT [pf] + dSe s [pS] + E v . (10) 

The first two contributions are the sums of the kinetic and strong-interaction 
energies and es[pf] is the surface energy density, which depends on all the 
particle densities in the surface layer _D S , and E v is the Coulomb interaction 
energy. The surface energy is given by integration of es over this narrow 
region around surface, where densities pf of, at least, some particle species 
change sharply. To be specific in this paper we concentrate on the discussion 



3 We consider here the case of the zero temperature and do not distinguish between the 
free energy and the energy thereby. 
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of the H-Q phase transition, as a concrete example. The corresponding values 
ejdn+str are introduced in Appendix A and the surface energy of the H-Q 
interface is discussed in Appendix B. 

We will further use that there is a shortest scale, d$, relevant for variation 
of above mentioned particle densities. In case of the H-Q phase transition the 
quark and nucleon fields change sharply and this scale is d$ ~ 1 fm, relating 
to the confinement radius (r\ ~ (0.2 0.4) fm ) and to the diffuseness 
of the nuclear layer (aV — 0.6 fm for an atomic nucleus). We will use 
that ds being much less than other scales involved in the problem, namely 
sizes of regions occupied by the phases and the Debye screening lengths 
Ag ~ (5 -T- 10) fm, see estimations below. Then one may replace the narrow 
surface layer by the sharp boundary 3D and approximate the surface energy, 
as we shall do, in terms of the surface tension a, taking J D dSes[pf] = 
J 9D dSe s [pf] = 4nR 2 cr for spherical droplet. In reality a is a complicated 
function of particle densities at the surface. It should disappear when energy 
densities of phases are equal, thus according to [27] producing a = as \ 
e 1 — e 11 | , with coefficient as slightly depending on particle densities of phases 
and e 1 ' 11 , as volume parts of energy densities of the phases. Simplifying, 
we will further neglect the density dependence of as(ps), considering as as 
constant. In this approximation there appears no contribution from the 
variation of a with respect to the pressure and to the chemical potentials, 
see (19), (21) below. The latter approximation is commonly adopted in 
literature, cf. [28] and refs therein. Since the values of as and a are rather 
poorly known we will allow their variation in the wide range. We will more 
closely discuss this question in Appendix B. 

The Coulomb interaction energy 4 Ey in (10) is expressed in terms of 
particle densities, 

1 f ^^QiPi^QiPj^) 



Ev 2 



(ii) 



with Qi being the particle charge (Q = — e < for the electron). 
Then the equations of motion (9) can be re-written as 

r)f a 

tf = ^§^-N?> a V a (r), N?' a = Q?/e, (12) 



4 We will further call it the "electric field energy" , in order not to mix it up with what 
we relate in this paper to "the Coulomb limit" with the spatially step-function dependence 

of pi. 
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with the electric potential well V a (r): 

J \ r-r> | " \ V ll (r), re D 11 [L6) 

generated by the particle distributions. 

We will keep the gauge invariance, so that V can be shifted by an arbitrary 
constant (V°) due to the gauge transformation, V(r) — > V(f) — V°. Formally 
varying eq. (12) with respect to V a (r) or pf we have the matrix form relation, 



do a 

A Z^k = N ?' a > ^B« k = 5 ik , (14) 



where matrices A and B are defined as 



X2rpa a a 

la _ -^kin+str r>a — U "i /-, r\ 

ij ~ 5pf5p° ' ij ~ dpf 1 ' 

Eqs. (14), (15) reproduce the gauge- invariance relation, 

clearly showing that a constant shift of the chemical potential is compensated 
by a gauge transformation of V a (r): fif — > p,f + AT ch ' a y°, as V(f) — > V(f) — 
V . Hence the chemical potential acquires physical meaning only after 
fixing of the gauge of V a {f); for the choice V° = (way I), /i" = /i" Gibbs , 
and for V° = -/i e , G ibbs (way II), \i« = 0. 5 

Applying Laplacian (A) to the l.h.s. of eq. (13) we recover the Poisson 
equation (re D a ), 

AV a (r) = 4ne 2 p ch ' a (f) = 4neQ?p?{r). (17) 

The charge density p ch ' Q (f) as a function of V a (r) is determined by the 
equations of motion (12). Thus eq. (17) is a nonlinear differential equation 
for V a (r). The boundary conditions are 

V l = V l \ W 1 = W 11 + 4vre 2 4 h , fedD, (18) 



3 The notation /if used in the paper [21] should be understood as M?Gibbs i n this sense. 
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where p§ h = a s h(S > ^ is the surface charge density. Below we will neglect 
a small contribution of the surface charge accumulated at the interface of 
the phases. Physically, typical scales for the spatial change of the charge are 
A^, 11 and R and, thereby, only a small charge can be accumulated at a much 
shorter scale d$, d$ < l/m n , = 140 MeV is the pion mass, the typical 
scale of the strong interaction. Variation of a over V yields a contribution to 
the surface charge Z s ~ Z dr ds/R at most that is <C Z dr , where Z dr is typical 
value for the charge accumulated in droplet, cf. eq. (Ill) of Appendix B 
below. 

We also impose the condition W u = at the boundary of the Wigner- 
Seitz cell, which implies that each cell must be charge-neutral. Once eqs. (17) 
are solved giving V a (r) and the potentials are matched at the boundary, we 
have density distributions of particles in the domain D a . 

Note that there are two conservation laws relevant in neutron star matter: 
baryon number and charge conservation. These quantities are well defined 
over the whole space, not restricted to each domain. Accordingly, the baryon 
and charge chemical potentials (//£ and /iq), being linear combinations of /if, 
become constants over the whole space, 



This fact requires two conditions for fif at the boundary dD, which determine 
the conversion of particle species of two phases at the interface. 6 

Once eq. (12) is satisfied, pressure becomes constant in each domain, 



-W = f Da df{e^ +str [pt (r)l - N?' a PmV a (r) - tfpftf} . (20) 



Hence, the condition of the minimum of f2 with respect to a modification of 
the boundary of arbitrary shape (under the total volume of the Wigner-Seitz 
cell being fixed) reads 



6 Each particle density is not necessarily continuous across the boundary, since it is only 
defined in each phase, while densities of leptons are well defined over the whole space. 
When particles of the same species i are allocated in both domains and the conversion 
of particle species becomes trivial, we must further impose the relations, fj,] = [i] 1 , and 
p\ = pY at the boundary. 



Pb = PB = Pq = Pq = HQ- 



(19) 



P l = 




(21) 
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where S is the area of the boundary dD, S = 4itR 2 , v l = 4irR 3 /3 for spherical 
droplet, and = in case of slab. Assuming as = const we dropped an 
extra small contribution to the pressure oc das/dps, see further discussion 
in Appendix B. 

The boundary of the cell does not contribute since all the densities are 
continuous quantities at this point. Eq. (21) is the pressure equilibrium 
condition between the two phases. 7 

The Debye screening parameter is determined by the Poisson equation, if 
one expands the charge density in 5V a {r) = V a {r) — around a reference 
value VJ.g f , which is also gauge dependent. Then eq. (17) renders 

ASV a (r) = 4ne 2 p ch > a (V a (r) = V T a e{ ) + (n a (V a (r) = V r a c{ )) 2 5V a (r) + (22) 

with the Debye screening parameter, 

= Cf)) 2 = 47re 2 

""(23) 

where we used eq. (16). Then we calculate contribution to the thermody- 
namic potential (effective energ y) of the cell up to 0(5V a (f)f. The "electric 
field energy" of the cell (11) can be written by way of the Poisson equation 
(22) as 

^ = / + / = / ^pgi dr - + I MI dr -, (24) 

Jd 1 Jd 11 Jd 1 87re 2 Jd u 8ire 2 

that is, in the case of unscreened distributions, usually called the Coulomb 
energy. Besides the terms given by (24), there are another contributions 
arising from effects associated with the inhomogeneity of the electric poten- 
tial profile, through implicit dependence of the particle densities on V I ' II (r). 
We will call them "correlation terms", u" OT = e£ in+str — /4*p"- Taking p" as 
function of V a (r) we expand e£ n+str in 5V a {r): 

4 m+s Mr-)]=e^ n+str (p t(K a ci )) (25) 

7 As was already noted, in a more detailed treatment of the problem with continuous 
density distributions, i.e. in the absence of a sharp boundary, the contribution of the 
surface energy is absorbed into P a . Hence P l = P u in such a more detailed treatment. 



dp' 



dV 



4 vrQf Q< 



dp? 



V a (f) = V a . 
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+ 



(//f + Nt' a V a (f 



1 

+ 2 



_a\2 



Aire 2 



+ (< + N?> a V"(f)) 



5V a {r) 

|2 



d 2 P? 
(dV a ) 



(5V a (r)) 2 + ... 



v«(f)=v£ 

We used eqs. (9), (13), (14) and (15) in this derivation. Using the expansion 

dpf 



-tifiW = -\rffi]vz t - 



dV a 



5V a (r] 



d 2 p a 



(dv a y 



(5V a (r)) 2 + ... 



(26) 



we obtain the corresponding correlation contribution to the thermodynamic 
potential fi cor = J D i druj l cov + f m dfuj^ OI : 

<r = ^i„ + str(P?(Cf)) - tfpUVS) - p ch ' a (Cf)Cf 



4vre 2 ' 8vre 2 ' '"' (2?) 

where we also used eqs. (22) and (23). In general V* e{ ^ V^ f and they may 
depend on the droplet size. Their proper choice should provide appropriate 
convergence of the above expansion in SV(r). Taking V^ ei = V^ f = V Te { = 
const we find 

(K a (Vr Ci )) 2 (V a (r)-V Ici ) 2 



8ne 2 



+ const, 



(28) 



and one may count the potential from the corresponding constant value. 
Here we also took into account that the term V r ^ { AV a (f)/(47re 2 ) does not 
contribute to Q according to the boundary conditions in the droplet center, at 
the droplet boundary (at zero surface charge), cf. (18), and at the boundary 
of the Wigner-Seitz cell. 



4 Conditions of chemical equilibrium. Charged 
particle densities for H-Q system 

Here we discuss the chemical equilibrium conditions on the example of the 
H-Q phase transition, although the scheme is quite general, as one will easily 
see. 
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The chemical potentials of charged particles (e and u, d, s quarks in our 
example) are determined by the equations of motion (9). Variation of (11) 
over pi gives extra 



SE l 



(29) 



contributions to the chemical potentials of charged particles. Using eq. (29) 
and eq. (93) of Appendix A, which determines the value eki n+str , we find 



Pq 



Pd 



= /i e =(3vr 2 pe) 1/3 + U I (f), 



1 + 



% a c\ 2/3 1/3 . l T/ I/^\ 



3tt J 



Ps 



1 + 



2a, 
3tt 

3tt 



m 1 It, * 



(30) 



The quark Fermi momenta are pp u 



Ps 



1 1/3 



(7T 2 p U ) 1/3 , PFd 



7r 2 p d ) 1/3 and 



PFs 



, a c is the QCD running coupling constant and m s is the mass of 



the strange quark. Note again that the chemical potentials are gauge variant 
before gauge fixing of V 1 (r). 

With the help of eqs (30), we further obtain 



Pu 

Pd 

Ps 
Pe 



1 



1 

1 - 



2a c 

7T 
7T 

2a ( 

7T 



2 l3 
Pu + ~V\r 



Pd ~ ^^(f) 



Ps 



-V\r) 
3 V ; 



3m 2 s p s 



= {p e -V\r)f/^\ 



(31) 



We use chemical equilibrium conditions for the reactions u+e <-> s, d <-> s, 
and n <-> p + e in each phase, 



Pu ~ Ps + Pe = 0, = p s , 



(32) 
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Pn pp Pet 



(33) 



and the conversion relation at the boundary, 



PB = Pn — 2pd + fi u , 



(34) 



which yield relations between quark and nucleon chemical potentials. 8 With 
the help of these conditions we obtain, cf. [17], 



Pu ^ 

Pd ^ 

Ps ^ 

Pe = 



1 

1 

7T* 



1 - 



1 - 



2a c 

7T 

2a c 

7T 

2a ( 

7T 



Hb + 2(-/i e + V\r)) 



Pb - {-Pe + V\r)) 



' Pb - (~Pe + V l (r)) \ 3 3m 2 s p s 



(p e - V\f)f/^\ 



(35) 



Since these formulae are invariant under the gauge transformation, V l (r) — > 
V l {r) — V° and accordingly fi e — > /x e — V by way of Eq. (16), we must fix 
the gauge to give a definite meaning of p e . In the way I discussed above 
in sect. 2 one puts V° = 0, and accordingly /i e = /i ej Gibbs- In the way II 
one takes V° = — /i e ,Gibbs, and p e = in this gauge. The expression for the 
electron density then renders p e = — (V I (r)) 3 /3ir 2 , cf. [18]. Thus, as we saw 
in sect. 2, V I (r) = — /4ioc(0 in this gauge. Note that these different gauge 
choices are reflected in the boundary conditions for V(r): V ll (r) — > V° as f 
approaches the boundary of the Wigner-Seitz cell, for the case f l — > (single 
droplet). 

For the quark matter with the number of flavors Nf = 3 that we study, 
we obtain from eqs (17), (31), cf. [17], 



Pch 



2a, 



TV 



X 



2fil(V\r) - fi £ 
9tt 2 



O 



'(v\?)-p e y 



Pb 



+ 



PBm; 
6n 2 



(36) 



8 Othcr conversion relation fi p — 2\i u + ^ is then automatically satisfied. 
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with the help of conditions (34). The electron contribution (oc (V l (r) — 
fi e ) 3 ) is numerically rather small compared to the quark contribution (for 
— V I (r) + /i e < m n of our interest) and can be neglected in the description of 
the quark phase. This means that the charge screening occurs mainly owing 
to the charged quarks rather than the electrons. Since —V l {r) + p e <C Pb 
we will also omit small (V l (r) — p e ) 2 contributions. Besides, m s ~ and 
the second term in squared brackets is essentially smaller than the first one. 
Numerically — V I (r) + p e ~ m n and the m 2 term gives a correction of the 
order of (V 1 (f) — p e ) 2 - Nevertheless, it is technically not difficult to keep 
it. Therefore, for generality we will retain below both mentioned terms, the 
linear in V l (f) — p e term and the constant one. 

The quark phases with Nf = 2 can be also realized in neutron star inte- 
riors in some models, cf. [22]. Therefore, we also present the charged quark 
density for Nf = 2: 



Here, the constant term is numerically not much smaller than the linear term 



For phase II, expanding the charge density Pch(r) around a reference value, 



* P*(V*if) = Kef) + *ft,(r) - Pe{V l \r) = Kef) - Spetf, (38) 



and using eqs. (14), (15), (23) and eq. (94) of Appendix A we find up to 
linear order 



faff * C Q \V l \r) - V£), Sp e (f) = ^-^(V l \r) - l£),(39) 

Co = ^, PFP = (^ Pp (v^f) = O) 1 / 3 , 



where A 2 2, and | A \ are the corresponding matrix element and the determi- 
nant of the matrix (15). 

Note that, at first glance, non-linear corrections to 5p p and 5p e are not 
too small, 8pf\r) = vr~ 2 (l / r I e I f — p e )(V u (r) — V^ f ) 2 and one may expect 




(37) 



in V 1 
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Sp^(r) ~ m n {V lI (r) — V^} { ) 2 . Nevertheless, these terms lead to rather small 
corrections of the resulting solutions V(r). The boundary conditions put ad- 
ditional limitations on values of corrections, that is often used in variational 
procedures and we have more, approximate solutions of the Poisson equa- 
tion. Finally, the correction to the resulting electric field profile due to the 
dropped nonlinear terms to Sp p (r) and Sp e (r) is proven to be on a ten percent 
level, cf. also [18], where the correction to Sp e (r) in an analogous screening 
problem was analytically evaluated resulting in only 1.5 % correction to V. 
We discuss the latter example in Appendix D. Explicit expression for the 
factor Co is derived in Appendix A. 

5 Spherical geometry 

5.1 Configurations of the electric field 

It is usually assumed that the thermodynamic potential should be invariant 
under the interchange of all the quantities related to the phases I <-> II, e.g. 
f l f U — 1 ~~ f 1 - Then one would deal with the quark droplets in the 
hadron surrounding, whether f l < 1/2, and with the hadron droplets in the 
quark surrounding, if/ 1 > 1/2. Even in case of unscreened configurations of a 
small size (in the Coulomb limit) this symmetry can hold only approximately. 
As we shall see, in the general case our equations are not invariant under 
these replacements and in order to conclude which structures are realized 
for intermediate values of f l in the vicinity of f 1 ~ 1/2 one should compare 
different variants. We postpone this study to the next paper. At least, 
for /' < 1 we deal with the spherical quark droplets embedded into the 
hadron surrounding and, at least, for 1 — f l <C 1, with the spherical hadron 
droplets embedded into the quark surrounding. Therefore, let us start with 
consideration of the spherical quark droplets of radius R immersed in the 
hadron surrounding. Increasing J 1 up to 1/2 we will check then the possibility 
of the possible transition from the structure of the quark droplets within the 
hadron phase to the structure of the quark slabs within the same hadron 
phase. 

The Poisson equation (22) with pl h from (36) describing the electric po- 
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tential of the quark droplet can be solved analytically. For r < R, we find 

V\r) - /i e = -^-sh(A) + Ul, (40) 

with an arbitrary constant Vq. For the Debye parameter k 1 and for the 
constant U\ we obtain: 

(^) 2 = ^ fl - — ) > (41) 



97T V 7T 7 ' u 4/lB 

For iV/ = 2 we would have 

Thus, the value Uq is rather small, especially for Nf = 3, and the main 
contribution to V l (f) comes from the first term in (40). Note that solution 
(40) is independent of the reference value V^ ci in this case, cf. (23), since p l ch 
in (36) is the linear function of V l (f) — /i e in the approximation used. 
Charge in the sphere of a radius r < R is given by 

Q(r) = rch(/tV) (\ - < 0, (43) 

being, thereby, negative, since Uq 1 > Uq and V l < 0. This negative charge 
is completely screened by the positive charge induced in the region R < r < 

For r > R, the Poisson equation with the boundary condition V'\r w = 
yields 

V l \r) -» e = V u *ch (k 1 \t - R w )) (1-5) + U U , (44) 
5 = th ( K n (%-r))/(K%), 

with an arbitrary constant V U , where the constant Uq 1 is given by 

U + fie = + Kef ( 45 ) 



22 



We take the reference value = V hnlk , where V hnlk is a constant bulk 
solution of the Poisson equation, (y hnlk — /i e ) 3 = — fxl Gihh8 = Sir 2 p p (V bulk — 
\i e = — // e ,Gibbs)> that coincides with the local charge- neutrality condition for 
the case of the spatially homogeneous matter. Note that V hnlk is a gauge 
dependent quantity, \/ bulk — ■> \/ bulk — V° at the transformation l /II (r) — > 
V n (r) - \/°: ^ = /i e , Gibbs , V bulk = for V° = and fi e = 0, V bulk = -/x e , Gibbs 
for V° = -/i e , G ibbs- 

The charge screening in the external region is determined by the Debye 
parameter 

{K y = 4e ^' - v "< » 2 + (46) 

7T C 

where the second term is the contribution of the proton screening. Taking 
p l % = 1.5po, /U e ,Gibbs — 170 MeV, /i B = /i n ~ 1020 MeV, a c ~ 0.4, we estimate 
typical Debye screening lengths as = 1/k 1 ~ ?>A/m n , and A" = l//t n ~ 
4.2/m T , whereas one would have A" ~ 8.5 /m n , if the proton contribution to 
the screening (46) was absent (Cq 1 = 0). With the estimate \% ~ A.2/m n 
we get that K n i? w > 1 for the droplets with the radii R > (2/ 1 ) 1 / 3 • 3.3/m^. 

5.1.1 Realistic case k R\\i > 1 

In this case eqs (40), (44) can be represented as 

V\r)-p e = ^sh(K l r) + U l , (47) 



V l \r) -p e = V^ch (k 1 \t - R w )) + (48) 
where Vq 1 ' 11 are arbitrary constants and 

i?w — -Rw ( I — rr-^ — arcth ^ 



Matching conditions yield 



^ (^-^(l + aoeth^Q) 
~ ch£ (a th£ • th(5i0 + 1) ' { } 
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y» „ (gj)Hgj (50) 

~ ch(5iO (a th£ • th(5iO + 1) ' 1 ; 

where we introduced notations 

«o = ^ = ^, Z = kR, ai = (/I)1/3 , (/)/ = — , (5i) 

and 

1 f(f) 1/3 \ 
5x = «i - -arcth — — . (52) 

? V «o4 / 

Tiny quark fraction volume. Let us consider the limiting case of a tiny 
quark fraction volume f l — > 0. Then, a given quark droplet does not feel 
the presence of other droplets which are placed at very large distances from 
it. The limiting expressions are recovered, if one takes i?w 00 • Then, 
solutions (40), (44) - (50) acquire simple explicit forms 

^(r) - ^ = t/ + (£/ 1 - £/) .; IB . im S 7 ^ (53) 

a sh{K L R) + ch^it) K x r 

'i _ thfcal fl 

y" (r) - „ e = £/» - ([/» - tfl X_ ^- exp (-«"(, - fl)) . (54) 

5.1.2 Less realistic case /t < 1 

In this case eqs (40), (44) can be represented as 

V\r)-^ = ^sh(K l r) + U l , (55) 

\/ n (r) - /i e = V^sh ( K n (r - i? w )) + (56) 
where V^ 1 ' 11 are again arbitrary constants and 

i? w = i? w fl tt— —arcth (k u R w )^\ . 
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Matching conditions yield 




(iff + ^11(5x0) 



(57) 



ch£ (a th£ • cth(a^) + 1) 




(gg - u£) (i - jthQ 



(58) 



sh(ai^) (aoth{ • cth(ai^) + 1) 



where now 



— -arcth 



( 



(/ I)l/3 



(59) 



At finite f l the Coulomb limit is recovered at ao£ <^ (Z 1 ) 1 ^ 3 with the help 
of equations (55) - (59), corresponding to the case « n i2 w < 1. However this 
limit is never realized for the realistic parameter choice, as we argue below. 

5.2 Effects of charge inhomogeneity in energy and ther- 
modynamic potential 

5.2.1 Realistic case k Ryj > 1 

"Electric field energy". We will calculate the contribution to the ther- 
modynamic potential (effective energy) of the Wigner-Seitz cell per droplet 
volume. Let us start with the proper "electric field energy" term ey = e v +€y: 




(1 + c^th^O) 2 (-|th 2 e + | - f th 2 £ + |thg) 
e (a th£ • th(3i0 + l) 2 



(60) 



being expressed in dimensionless units (51) and 




(61) 
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where we used eqs. (24), (40), (49). With the help of eqs. (44), (50), from 
eq. (24) we find 



v 4itR 3 Jr 87re 2 



^ (l ~ |tK) 2 (l - |aoaig(l ~ thVQ) + looethgiQ) 
~ e («othC • th(5i0 + l) 2 

a 2 (l-ith^) 2 (l-th 2 (5i0) 
_/3 °2(/i)V3 (a§f(/V /s - 1) («oth£ • th(5i0 + l) 2 ' (62) 

As above, limiting expressions for the case of a tiny quark fraction volume 
f 1 — > are reproduced with the help of the replacement i?w - ► oo. 

"Correlation terms" to the thermodynamic potential. Recalling that 
the solution V 1 is almost independent of the reference value V^ ef , we choose 
Kef = Kef = V hnlk in order to explicitly calculate correlation terms; then 
8V l (r) = V^r) — V bnlk for r < R. Averaging (28) over the droplet volume, 
with the help of (40), (49), we obtain 

c i _ 3 f\^,.i „ & (l + ^th^)) 2 



- ----- f 47rrW ~ t h£ - - , 

2/3 (l + aoethCaiQ) ft 
+ ^(a the-th(5 1 + l) ( H ^ ~3 ' ( j 

In the hadron phase introducing SV u (r) = V ll (r) — \/ bulk , and using eqs. (28), 
(44), (46), (50), we find 

2)1 = — ^ / Anr 2 drul r ~ (64) 

./ft 



, 2 



(l - |th£) [|th(5x0 + c*i(l - th 2 ^)) 



2 (a th£ • th(5i0 + l) 2 

Poa 2 (l-lth^'q-th 2 ^)) 
2 (f) i/3 (« th£ • th(5x0 + I) 2 KW)- 2/3 - 1) ' 



Without the (^-correction term in (44) we would come to the very same 
equations but with 5i replaced by «i and without the last terms in expres- 
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sions (62), (64), cf. [21]. 9 The latter terms disappear in the limit f l — > 0, 
a.\ — > ai. Thus, account of the 5-correction results in a small correction, 
at least, in both limiting cases cto^/ 1 ) -1 / 3 3> 1 (small quark fraction) and 

'>,,((/'! ; :i L« I- 

We could also use other values for V^, f and V^ f , e.g. we have checked that 
using general eq. (27) with V^ c{ = V^O) and V^ f = V u (Rw) leads to the very 
same result. These two choices are very natural although one could also select 
V^gf 1 differently varying their values within interval (\/ bulk , 0). Corrections are 
expected to be of the same order of magnitude, as those coming from the 
above dropped non-linear terms in particle densities. 

5.2.2 Less realistic case, k -Rw *^ 1 

In this case eqs (60), (62), (63), (64) are rewritten as 

6y ~ 4ttWo 87re 2 47rr dr 

(1 + ao ecth(5 1 Q) 2 (~|th 2 e + | - |th 2 C + |thj) 
e (a th£ • cth(5i0 + l) 2 ' 



v 4ttR 3 Jr 8ne 2 

(l - |th^) 2 (l - |a ai£ 2 (l - cth 2 (5x0) + ±a fcth(aiO) 



- A 



0" 



£ 2 (a th£ • cth(5i0 + l) z 



qg(l-lth^) (l-cth'fofl) 

Mo 2(/i)V3 (i _ c^ea 1 )- 273 ) (« tK • cth(5i0 + 1) 2 ' 

for the "electric field energy" terms and 

~J - /VrW ~ Ac (l + ^cth( 5l Q) 2 / _ J_\ 

2/3 (1 + Q!o^cth(aiO) , . & 

+ e (aoth^.cth(5iO + l) W j 3' 10 



9 Please pay attention on two misprints in [21]. In eq. (36) factor 3 is missed and in eq. 
(39) the term 4/3 should be replaced to 1 
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^cor 



Poa i 1 - | th £ 



f Rw 1 TT 

J Anr druj COT ~ 



(68) 



kth(5i0 + cti(l - cth 2 («iO) 



2 



(a th£ • cth(5i0 + 1) 

l-ithe) 2 (l-cth 2 (5i0) 



2 (/i)i/s ( ao the • cth(5 1 + l) 2 (1 - aieU 1 )- 2 '*) ' 
for the correlation terms to the thermodynamic potential. 

5.2.3 Surface energy 

In our dimensionless units the total quark plus hadron surface contribution 
to the energy per droplet volume renders 



es/A) = Pi/i, Pi = 3kV/A), 



(69) 



see (61), and we used that es = 3a/ R, cf. eq. (110) of Appendix B. The 
coefficients p , Pi are evaluated with the help of eqs. (41), (61) and (69). For 
N f = 3, 



(70) 



A^(^-^(l-^)" /2 ^^^^^. (71) 

For the typical values -£/ n ~ /^.Gibbs ^ 170 MeV, /x B ~ 1020 MeV, a c ~ 0.4 
and m s ~ 120 -j- 150 MeV we estimate po ~ 1.6m£. Thus, with the value 
a ~ 1.3m^ we obtain Pi ~ 0.7, whereas with er ~ 10 MeV/fm 2 ~ 0.14m^ 
we would get p\ ~ 0.08. 

Note that all the energy terms at £ = £ m in, i.e. at the energy minimum 
point, are expressed in terms of only four parameters: «o = Afj/Af], f l -> Po 
and Pi, and the dependence on the properties of the concrete system under 
consideration is hidden only in the values of these parameters. 
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Correlation term to the surface energy. In the general case, there is 
still a correlation term to the surface energy due to a dependence of a(rii(V)). 
However this term is rather small being proportional to the variation of the 
surface charge SZg/SV. The physical reason of this smallness is obvious. 
Electric fields are changed on large scales (A^, 11 , R) and cannot substantially 
affect surface quantities which are changed at shorter distances (~ <ig). We 
have omitted a contribution of the surface charge to the charge distributions 
and by the same reason we omit the correlation term to the surface energy. 

5.2.4 Coulomb limit 

Tiny quark fraction volume. Peculiarities of the Coulomb limit at finite 
values of the quark concentration f 1 are discussed in Appendix E. For the 
case of a tiny quark fraction volume (very small f 1 ), i?w is very large. Then 
one puts pl h = const, p 1 ^ = within the Coulomb limit of Appendix E and 
recovers the corresponding expressions. 

/,From our general equations describing inhomogeneous charge profiles we 
reproduce the Coulomb limit for the case of a tiny quark fraction volume, 
if we first put a\ — > oo, and then expand the terms e\ + ty + eg in £ <C 1. 
Thus, we should put £ <C 1, «i£ ^> 1. Then we, indeed, recover the Coulomb 
plus surface energy per droplet volume 

?c,s = ?c + e s = ^ (^£ 2 + i£ 2 + |), (72) 

where partial contributions correspond to the terms e v , e v l and eg. We 
needed the Taylor expansion of functions entering e", see (62), up to two 
non- vanishing terms in order to recover the contribution oc £ 2 and we needed 
the Taylor expansion of functions entering e v , see (60), in £ 1 up to three 
non- vanishing terms. 

Both the correlation terms uo l COT oc £ 4 and u)" r oc £ 3 can be dropped in the 
Coulomb limit, for droplets of a tiny size £ «C 1. For u)* or , see (63), the Taylor 
expansion in small £ should include four non-vanishing terms to recover the 
term oc £ 4 , and one needs to do the Taylor expansion of the functions entering 
(64) for ojUj. up to two non-vanishing terms with the result oc £ 3 . 

With the help of eqs (36), (41) the value e c in (72) is easily recasted in 
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the standard form, 

ec = ^f{ P \?R 2 - (73) 

The function ec,s has the minimum at £ = £ m in, corresponding to the optimal 
size of the unscreened droplet: 

e C) s(£min) = 36 C (Un) = ^ , £mi„ = ' ■ (74) 

The Coulomb limit is recovered only for £ min <C 1, i.e. p 1 <C 0.1, whereas 
with above estimate (3i > 0.1, cf. (71), we always get £ min > 1. Thus, we 
conclude that the pure Coulomb limit is never realized within mixed phase. On 
the other hand, the criterion of the tiny quark fraction volume ai£ m in 3> 1, 
which we have used, is applicable only for f 1 <C chqPi. 

Finite quark fraction volume (f 1 ) in the Coulomb limit. As above, 

we assume £ min ~ p\ /S < 1, but now a fmin < (f) 1/s (0.5 > f 1 > ag/3i). 
Then, from our general equations (65), (66) we find 

y + ey " is^ + ^i-/ 1 )) 2 ' (75) 

that coincides with the asymptotic form expression (142) of Appendix E for 
the Coulomb energy obtained there with the step-function charge densities. 
Again this limiting case is not achieved for realistic parameter choices, since 
then we always have £ m i n £ 1 and ao£min ~ 1- 

5.2.5 Limit of extended quark droplets 

In the limit «i£ ^> 1, £ 1, corresponding to the large size droplets, from 
(60), (62), (63) and (64) we find that all the terms contribute to the effective 
surface energy density (neglecting the curvature terms oc l/£ 2 ), 

~I _ A)«Q 1 ~n A)<*o 1 

e V 0/1 , „ \2 <■■> e V 



2(l + a ) 2 e' V 2(l + a )H' 



W cor 



2{3 a /3o«o 



l + ao) 2(1 + a o y 



£ + T' c ° r 2(l + a )H' ( } 



30 



Thus, 5u tot /f3 = (e y + £* or + + e s )/(3 = Pf^ 1 + const with pf = 
S^alrt* 1 1 ' (5q and therefore, in the given limiting case of extended droplets 
the electric field effects can be treated with the help of an effective surface 
tension. The full surface tension o"to t hcr then renders 

sphcr „ \I Po a o[ a + !] 

^tot = <J + a v = a - A D — . (77) 



3(1 + a ) 

The first a term is the contribution (related to the scale ds) of the strong 
interaction. The second negative term is the contribution of electric field 
effects (related to the scales A^ n ) . It depends largely on the values of param- 
eters. For /i e , G ibbs = 170 MeV, \i B = 1020 MeV, m s = 150 MeV, a c = 0.4, 
we estimate the contribution to the surface tension from the electric effects 
as a v ~ -60 MeV/fm 2 . Thus, for a < a c ~ 60 MeV/fm 2 we deal with the 
mixed phase, whereas for o > a c , with the Maxwell construction. Note that 
only due to the correlation term u) l COT we obtained negative contribution to 
ay. Also note that the spatial variation of the electric field within the surface 
layer of thickness ds may only slightly affect the "old" value a at ds A^ 11 - 
The dropped ~ l/£ 2 terms are important to determine the size (£ m in) of 
the droplet within the mixed phase. For very large droplet radius these terms 
can't compete anymore with ~ l/£ terms. For such large radii of droplets 
the mixed phase already can't exist. With this limiting case we may describe 
the charge distribution within the Maxwell construction. 



5.3 Numerical results 

In Fig. 2 we demonstrate the dependence of the contribution of the inho- 
mogeneous charge distributions to the total thermodynamic potential per 
droplet volume for the case of the lattice of spherical droplets, 5u to t/Po = 
(ey + €y + u l cor + ul l or + ?s)/ A), given by the sum of partial contributions (60) , 
(62), (63), (64) and (69), if k, u R w > 1, and by the sum of (65), (66), (67), 
(68) and (69), if k u Rw < 1. It is presented as the function of the droplet 
size (£ in dimensionless units) for two values of the Wigner-Seitz parameter 
f 1 = / = 0.01 (dashed curves) and / = 0.5 (solid curves) at fixed value 
Pi (each panel) for a = 1, as the representative example. The curves for 
/ = 0.01 deviate only little from those for / — > related to the single droplet. 

The curves labeled by "C" , the Coulomb curves, demonstrate the contri- 
bution in the Coulomb limit: ec,s/Po is determined by eq. (72) for / — > 
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and by eqs. (142), (147) in general case, dashed curves stand for / = 0.01 
and solid ones, for / = 0.5. Each "C" curve has a pronounced minimum 
a t £ — £min oc P^ 3 . For £ > £ min the "C" curve shows quadratic growth 
deviating drastically for £ > 1 from the other curves. 

The dashed and solid "e.m." curves shown in each panel demonstrate the 
contribution e e . m .//3o = (ey + £s)/Po f° r / — 0.01 and / = 0.5 respectively, 
ignoring correlation terms. The quantity e e m . for / — > is the counterpart of 
that in the Coulomb limit of ?c,s, when the charge screening effect is taken 
into account. We see that the minima at the "e.m." curves disappear already 
at Pi > 0.03. 

For Pi < 0.01, minimum points of the dashed "C" curves £ m i n deviate 
little from the minima of the dashed "e.m." curves. Only for such small 
values of Pi and £ min we recover the Coulomb limit! However, one may 
obtain such small values of Pi only for tiny values of the surface tension 
and very large values of the baryon chemical potential. With increase of the 
latter, the Debye screening parameter k 1 is also increased and the droplet 
radius R m i n = ^min/^ 1 is proved to be essentially smaller than l/m n . For 
larger values of Pi, the difference between the minima of the dashed "e.m." 
and the dashed "C" curves is proved to be more pronounced. The difference 
between the minima at the dash "C" and dashed curves ( 5uj t ot/Po) for the 
case of a tiny quark fraction volume (/ = 0.01) is minor for p 1 < 0.05 
but becomes pronounced for larger values Pi, whereas for the solid curves 
(/ = 0.5) this difference is always pronounced, even for very small Pi. The 
latter is due to an essential contribution of the correlation effective energy 
in this case (/ = 0.5), whereas p\ h is fitted by coincidence of the "C" and 
its generic partner the "e.m." curve in the limit ^<C 1. In the framework of 
the Coulomb limit the value p* h is defined in the variational procedure: the 
total energy is assumed to be dependent on p l ch via / and p\ h is determined 
by the minimization of the energy. Then mentioned discrepancy could be 
partially diminished. However even by normalizing the solid "C" curve to 
coincide with the solid curve at £ 1 (for that please simply do the shift of 
the minimum at the solid "C" curve up to the value given by the minimum 
at the solid curve in the panel Pi = 0.01 and use the same value of the 
shift considering other panels) we obtain that the difference between the 
minima at the curves (which determine droplet radii) becomes pronounced 
for Pi > 0.05. This difference is more pronounced for a 7^ 1- 

The dashed and solid curves converge to 1/3 with the increase of £. The 
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large £ asymptotic of the curves (counted from 1/3) is oc l/£, being inter- 
preted as the surface energy term, characterized by a significantly smaller 
value of the surface tension (77) than that determined only by the strong 
interaction. We see that for a >| oy |, corresponding to (3\ > (5 lc ~ 0.5 for 
ao = 1, see Fig. 3, the structured mixed phase is proved to be prohibited, 
since a necessary condition of its existence (the presence of a minimum in 
the droplet size) is not satisfied. Contrary to this in the Coulomb limit the 
necessary condition is always fulfilled. Large-size droplets are realized within 
the mixed phase if a <\ ay |, | a + ay |<c| ay | and within the Maxwell 
construction, if a >| ay \. 

The difference between the solid and dashed curves and the corresponding 
"e.m" curves shows the important contribution of the correlation energy in 
the H-Q structured mixed phase. This correlation contribution contains two 
parts, the model independent constant term (/3o/3 in our case), as follows 
from (63), and a model dependent part. The latter however depends only on 
three parameters /?o/3 and «o — ^d> ^d> The constant term /?o/3 shows a 
difference in the energies of the Maxwell construction and the mixed phase, 
if one neglected finite size effects in the calculation of the latter. 

The minima of the solid "e.m." curves are always below the corresponding 
minima of the dashed curves, however for the 5u)tot(£min)/A) at the dashed 
and solid curves the situation changes to the opposite one, &Dtot(£min)/A) 
increases with /. Comparison of the solid and the dashed curves for fa > 0.5 
demonstrates only a minor dependence of the value 5uj tot //3 on the Wigner- 
Seitz parameter / in this case. 

Fig. 3 demonstrates the dependence of the droplet radius £ m ; n (in di- 
mensionless units) on the value f3\ for different ratios of the Debye screening 
lengths ao, again for two values of the Wigner-Seitz parameter / = 0.01 (a 
tiny quark concentration ) and / = 0.5. As seen from the figure, the depen- 
dence of £ m i n on the parameter a is rather sharp, whereas it was completely 
absent for the Coulomb curves in previous works since there p l ch was consid- 
ered as a function of / rather than of k 1 and k 11 . The dependence on / is 
more sharp for a < 1 and less sharp for cto > 1. All the curves demonstrate 
the presence of a critical parameter set (/3i c ,£mm)- The values f3i c are shown 
by arrows in Fig. 3, cf. with Fig. 2. Corresponding values £^ in are rather 
large demonstrating that even the droplets of the size > 10A}-, may still exist 
within the mixed phase. This statement disagrees with what was expected in 
previous works using the Coulomb solutions. It was thought that the mixed 
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phase should disappear for £ > Ad due to the screening. 

For larger values of (3\ {(5\ > f3i c ) the Maxwell construction becomes 
energetically favorable. The Coulomb curves, "C" show the value £ min related 
to the unscreened the function of U (Um(Pi) = (15/V4) 1 / 3 for 

a tiny quark fraction volume). The necessary condition for the existence 
of the mixed phase (existence of £ m i n ) is always fulfilled in the Coulomb 
limit, opposite to what we find in our general consideration. Comparison of 
the curves with the "C" curves again demonstrates the importance of the 
screening effects. 

6 Plane geometry 

6.1 Electric field configurations in Wigner-Seitz cell. 

Without incorporation of screening effects previous papers argued that with 
increasing of the quark fraction volume the lattice of spherical droplets first 
undergoes the phase transition transforming to the structure of rods and 
then to the structure of slabs, which eventually become energetically more 
favorable compared to spherical droplets and rods at f 1 < 0.5, cf. [16, 17]. 
Therefore, we will describe these configurations but now including the charge 
screening effects. We will concentrate on discussion of slabs, since in absence 
of screening effects, as has been shown in mentioned works, these structures 
appear already at sufficiently small quark fraction volume, and also, since in 
this case we may easily proceed further with our analytical approach. 

Let the quark slab occupies a layer x G (—R,R), — oo < y,z < +oo, 
being placed in the Wigner-Seitz cell x G (— Rw, Rw), Rw > R- We suppose 
that in the interval x G (-R, R) the electric potential well V is an even 
function of x. Also we assume that the solutions of eq. (17) should not 
change under simultaneous replacements x — > — x, R — > — R, -R\v — — -Rw- 
For the charged density (36), dropping the terms ~ 0(V 2 (x)), we easily find 
a solution of (17) 

V\x) -n e = V$ch(K l x) + U\, —R < x < R. (78) 

Here, by superscript (1) we indicate the flat geometry of the structure, Vqi is 
an arbitrary constant and the values k 1 , and Uq are determined by eq. (41), 
as in the case of the spherical geometry. 
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The solution of (17), V u {x) = \/ Q n + 5V u {x), Vj 1 = const, | 5V u {x) |«| 
Vq 1 I for R < x < i?W) with the boundary condition V' x {x = Rw) = 0, reads 

V n (x) -n e = V$ch(K U (x- R w )) + Ul\ R<x<R w , (79) 

where Vq n is an arbitrary constant and the values n 11 and Uq 1 are determined 
by the same equations, as in the case of the spherical geometry. The corre- 
sponding solution (17) for — Rw < x < —R, with the boundary condition 
VJyX = —Rw) = renders 

V l \x) - /i e = V$ lC h(K n (x + R w )) + U$, -R w <x<-R. (80) 

Matching of potentials and their derivatives at x — ±R yields 

yd) = (ffi ~ U>0 (gl) 

0,1 oochZ + sh£cth(o: 2 £) ' 



yd) = (ui - u$m (82) 

0,11 a sh(a 2 + th£ch(a 2 ' 1 ' 

where for one dimensional structures the parameter <y,\ is replaced by a 2 = 
a (l-f)/f\f = R/R w . 

Boundary conditions for x = —R, x = —Rw are automatically fulfilled. 
The solution describing a tiny quark fraction volume, when one may deal 
with a single slab, is obtained from here in the limit Rw —> oo. 

6.2 Effects of charge inhomogeneity in energy and ther- 
modynamic potential 

Using (78), (81), we first find contribution to the energy per slab from the 
term eyj: 

, a > = _l [> <?xML ix = A ^(^)K-i + |ty) (g3) 



2R J-R 8ire 2 



6(a 1 th^ + th(a 2 O)' 
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Using eqs (79), (82), we then obtain the term eyjj: 



(1) _ 1 (VVW 



B th 2 e(^th(a 2 Q + g(th 2 (a 2 0-l)) 
6 (c*o Hh£ + th(a 2 0) 

For £5^, with the help of eqs (28), (78) and (41) taking = U bulk we find 

1 ,R th 2 (q 2 Q (l-th 2 e + |thC) 

2K 6 («o + th(a 2 0) 

2th(a a flth£ ft 

Using eqs (28), (79), (46) and taking = U bulk , we obtain 

th 2 e [|th(a 2 + a 2 (l - th 2 (a 2 0)" 

= ^° r - 1 \2 — ■ 

6a [ a o th£ + th(a 2 e)J 



The surface energy contribution per slab is given by 



see eq. (110) of Appendix B. 



4" = (87) 



Tiny quark fraction volume. Single slab. Expressions for a single slab 
are obtained taking the limit i?w —> oo (a 2 — > oo). 

Slab of a tiny transverse size. The solution for a single slab of a tiny 
transverse size is recovered for £ <C 1, a 2 £ ^ 1- This limit is never realized for 
an unscreened slab, the potential of which is linearly growing for x — > ±oo. In 
the realistic case under consideration any charged slab is screened already at 
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distances ~ from it. Then, the main contribution to the energy per slab 

comes from the terms eyjj ~ cDcor.n °^ £■ The term lyj is oc £ 2 , and £>cor,i oc £ 3 . 
In this limit, i.e. for a tiny quark fraction volume, we immediately obtain 
that unscreened spherical droplets are energetically more favorable compared 
to slabs. Indeed, we have 

eg> s ^ M + (88) 
c ' s 3a 3f ' K ' 

from where we find the minimum £ min = and ec,s(£min) = 2ec(£ m in) = 

Pi/ot-o, that is larger than the term (74) for the case of the spherical 

droplet (ec,s(£min) oc [3^ 3 ) for f3\ < 0.11/ccq. This statement holds in the limit 
f3ia <^ 1 and f 1 a (/3iao) 1 ^ 2 corresponding to £ min 1 and a^min ~> 1 
for slabs, and fa <C 1, af,Pi ^> f 1 corresponding to £ min <C 1 and ai£ m i n 1 
for droplets, accordingly. 

Slab in Wigner-Seitz cell, both of tiny transverse sizes. This limit 

is recovered for £ < 1, a 2 £ 1. Then we arrive at expressions 

~<i) _ alalM 2 ~<i) _ ayOtlM 2 , RQ x 
^-9(1 + 00^)=" ^ n -9(l + a « 2 ) 2 ' ( j 



.(i) _ A) 2(3 a al(a + q 2 ) 

;cor ' 1 ~ 3(1 + a a 2 ) 2 9(1 + a a 2 ) 3 ' 



„(i) Poaoa 2 2p a al(a + a 2 ) 2 / n nA 

3(l + a tt2) 9(l + a «2) 

and thus oo^ Tl + u;£2,ii ~ §A)(1 + ckoC^) -1 . 

Summation of all the terms (89), (90) together with (87) yields, after 
minimization over £: 

uV>(P )/3 - £ /3 Ml-f) 1/3 | f m) 

^totUminj/A) 3 l/3 2 2/3[ a 2 + _ ^2/3 + 3^2 + _ ^2)] > l yi J 

where we used that 

Un " I 2ag(l - f) ) <<h (92) 

and a2^mm <C 1, that corresponds to (/^ao) 1 / 3 ^ 1 and J 1 ^> ao^i^o) 1 / 3 . 
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Single slab of large transverse size. As follows from eqs (83) - (87), 
in the limit of a single large slab of a much larger transverse size than the 
typical Debye screening lengths (£>■!, a 2 £ 3> 1), all the terms contribute 
to the surface energy per slab with the corresponding surface tension terms 10 
= o"| phcr , and, thus, for the total surface tension we get a[H = o"tot her - 

The slab of a large transverse size has three times smaller negative contri- 
bution to the surface energy per slab than that for the corresponding spheri- 
cal droplet in case Otot < 0> when the mixed phase is energetically favorable. 
Here, the factor three comes from the space dimension. Thus, droplets are 
energetically favorable compared to slabs also in this limit. 

For the case a[H > realized within the Maxwell construction the plane 
boundary layer is energetically preferable. In spite of that, in any case the 
boundary layer, which arises between two separate phases within the Maxwell 
construction, is spherical, following the geometry of the neutron star as a 
whole. Above discussed slabs could exist in the neutron star only, as config- 
urations of the mixed phase, being rather small in their transverse sizes. 



6.3 Numerical results 

In Fig. 4 we show the same quantity, 5uj to t/(3o, as in Fig. 2, but now 
for slabs, for the values of the Wigner-Seitz parameter / = f 1 — 0.01 and 
/ = 0.5, and for the parameter a = 1- The "C" curves are calculated with 
the help of eq. (88) for / = 0.01 (dashed) 11 and eqs (143), (147) of Appendix 
E for / = 0.5 (solid). Eq. (143) gives a divergent result in the limit / — > 
and is not applicable to describe single slab of a tiny transverse size. 

Difference between the non-labeled curves (Su)^/ /3q) and the "e.m." curves 
shows how large are correlation terms. We see that forms of all the curves 
are analogous to those for the spherical droplet case. Minima disappear at 
the "e.m." curves already for f3\ > 0.03 at / = 0.01 and for f3\ > 0.05 at 
/ = 0.5. Minima of the non-labeled curves (Sul^; / Pq) disappear at (3\ > 0.5 
for both / = 0.01 and / = 0.5 cases, cf. Fig. 5. 

10 The curvature terms oc l/£ 2 are different in case of droplets and slabs. 

11 Thus, this is actually not the Coulomb curve, being calculated without inclusion of 
screening, rather it is the "C" curve obtained from our general equations in the corre- 
sponding limit of a tiny quark concentration and £ <C 1. We call it "C" in the sense that 
expression (88) is applicable for £ <C 1 but we apply it at arbitrary £, as it is done with 
the Coulomb solutions. 
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The dashed " C" curve is the appropriate limit of the dashed curve (Sutot/ A)) 
at a tiny / and 1. It grows linearly with £. At their minima the dashed 
" C" curves begin to essentially deviate from the minima of the corresponding 
dashed curves for (3\ <; 0.05 -j- 0.1. If they were normalized to coincide with 
the "e.m." dashed curves at a tiny /3i, their minima would be shifted down 
by factor l/y/2 since ey n ~ w^iii ^ an d ^vi ^ £ 2 an d ^iori ^ £ 3 are 
neglected. Solid "C" curves (/ = 0.5) are normalized to coincide with the 
corresponding "e.m." curve for £ <§C 1, see discussion in Appendix E. 

Fig. 5 demonstrates the dependence of the slab size in dimensionless 
units, £min, on the parameter f3\ for / = 0.01 and / = 0.5, for three values 
of cto, cf. with Fig. 2 for the spherical droplet case. The ao dependence is 
here even stronger than for the spherical geometry. Also the transverse size 
of the slab is smaller. 

7 Comparison of energies of droplets and slabs 

As has been argued in the previous studies, cf. [8] sect. 4.3.1, the bulk energy 
densities of the quark and hadron phases are about two orders of magnitude 
larger than the sum of the Coulomb and surface energies. This large factor 
is partially due to the fact that the latter sum contains an extra oc e 2 / 3 pre- 
factor compared to the energy densities of the bulk phases given by eqs (93) 
and (94) of Appendix A. One can see this pre-factor after minimization of 
the sum of the surface and the Coulomb energy densities with respect to R, 
taking into account that the values p^ 1 entering eq. (140) of Appendix E 
do not contain this pre-factor and that the resulting surface energy density 
is of the order of the Coulomb one. An extra small factor arises from the 
smallness of the difference p* h — p l ^ h entering eq. (140). 

With our screened potentials, the smallness of the value 5a3 to t(£min) is due 
to the parameter A oc e (cf. eq. (71)) for typical values £ m i n Z 1. 

/,From the Figs 1-5 and eqs. (93) and (94) of Appendix A we may see 
that typically 5cDtot(£min) (with £ min essentially dependent on the geometry 
of the structure) is at least by an order of magnitude smaller than the bulk 
values of the energy densities. This justifies the same "two part approach" to 
the problem, as has been used in previous papers of Glendenning and his co- 
authors, of computing first the bulk properties and then on this background, 
the geometrical structure imposed by the effective energy &^tot(£min)- 

The bulk structure is defined following the strategy, sketched in [8], sect. 
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4.1.3. Begin with solving the equations that define quantities of the phase II 
using the local charge neutrality condition and increasing the baryon density. 
At each baryon density find the values of the chemical potentials /x§ )loc and 
/ig I loc . At each density also find the solution of the high-density phase I at the 
very same value of the chemical potentials thus putting p}^ loc = pB,Gihbs — 
const and p^\ oc = /i e ,Gibbs = const and using the relations (32) - (34). Locate 
the value of the chemical potentials for which the pressures in the two phases 
are equal. The beginning point of the mixed phase is fixed in this way. The 
complementary procedure will locate the boundary between the mixed phase 
and the pure phase I (quark phase). Then we may find the properties of the 
mixed phase in which both phases are present in equilibrium. To find the 
volume fraction of phase II at fixed ps, for given / (starting from infinitely 
small /) solve the equations defining both phases subject to the conditions of 
equal pressure and global charge neutrality, where the charged densities are 
defined through the values of the corresponding constant chemical potentials 
which, as well as the pressure, vary with /. The value / is a function of 
Pb, see (5). The increase of the value / corresponds to the increase of the 
baryon density (5). The value / = 1 corresponds to the end point of the 
mixed phase. 

Then, in order to define, which geometrical structure is energetically fa- 
vorable, we may use the values of ps and / thus calculated, and compare the 
quantities ^tot^(^n) an d &^tot v£rrdn)- The latter values relate to the dis- 
tinct possible configurations (quark/hadron) of different geometries d\ and 
d 2 . As we have mentioned, in the present paper we restrict ourselves to 
the consideration of the quark structures embedded into hadron matter for 
/ < 0.5 and by studying of the two geometries di — 3 (spherical quark 
droplets) and d 2 — 1 (quark slabs). Each quantity &*)tot (£min) and <5^tot(Cnin) 
is minimized with respect to the structure size (£^ n and £mhJ- The values 
and are, thereby, essentially different, see Figs 3 and 5. This is due 
to the fact that the minimization in £ is performed for the quantities Suj^ t (£) 
and &t)tot(£)> bd n g of the same order of magnitude (having the same e pre- 
factors). Thus the argument which allowed us to fix / to be the same for 
both structure geometries does not hold for £. As we have demonstrated 
above, cf. Figs 3, 5, the minima and exist not for all values of 
Pi, but only for those related to not too large values of the surface tension 
parameter. Otherwise the presence of the mixed phase is not energetically 
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favorable compared to the configuration given by the Maxwell construction. 

Above procedure implies that the corrections to / due to its dependence 
on the geometrical structure and accordingly on the values of £^; n or £® n 
are rather small. Therefore they are neglected. The values of p^ 1 used in 
determining the parameters of the inhomogeneous spatial structures, should 
not be constants in the general case, rather they are now functions of r given 
by (36) and (38). The global charge neutrality condition is now satisfied 
since we use the boundary condition VV = at the boundary of the Wigner 
- Seitz cell. 

Fig. 6 shows which structures (spheres or slabs) are more energetically 
favorable for the given value of the parameter /, in dependence on the value 
of Pi. The solid curves, Su^/Po = &*5tot(£min)/A), present droplets and 
the dashed curves, 5ui™ t n /p = S&ill (£min) / Po, the slabs. The result, which 
structure is energetically favorable, essentially depends on the values of pa- 
rameters cto an d Pi- For a — 0.5, Pi > 0.3 spherical droplets are energetically 
favored for any value of the parameter / < 0.5. For a — 2, Pi Z 0.3 slabs 
are energetically favored. For unrealistically large values of ceo, see the cases 
a = 5; 10 in Fig. 7, we find that slabs become to be energetically favorable 
structures already for Pi > 0.01. It occurs then for any /. Our conclusions 
here essentially differ from those derived in the Coulomb limit, if one ignored 
the screening effects for any Pi and assumed the step-function charge- density 
profiles. In the latter case the droplet phase is proven to be energetically fa- 
vorable for / < f v c oA ~ 0.2, as has been shown in ref. [16, 29]. Then for 
fc° d < f < fc lab — 0-32 rods are prefered and for / > f* lab slabs are prefered. 
We recover this limit case result from our general expressions as well, but 
only for tiny values of Pi and = 1, see discussion in Appendix E and 
Fig. 8. In the latter we present the effective energies of droplets and slabs 
obtained with our general expressions in case Pi = 0.001 for several values 

of CKq. 



8 Concluding remarks 

Summarizing, in this paper we have discussed the possibility of the pres- 
ence of structured mixed phases at first order phase transitions in multi- 
component systems of charged particles. Using the thermodynamical poten- 
tial as the generating functional we have developed a formalism, which allows 
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to explicitly treat electric field effects including the Debye screening. 

We have studied a "contradiction" between the Gibbs conditions and 
the Maxwell construction extensively discussed in previous works. We have 
demonstrated that this contradiction is resolved if one takes into account 
the difference in the treatment of the chemical potentials used in the two 
approaches. The different values of the electron chemical potentials in the 
Maxwell construction and the ones used in the Gibbs conditions do not con- 
tradict each other if one takes into account the electric field arising at the 
boundaries of the finite size structures of the mixed phase and at the bound- 
ary between phases within the Maxwell construction. 

As an example, our formalism was applied to the hadron-quark struc- 
tured mixed phase. Using a linearization procedure we analytically solved 
the Poisson equation for the electric potential and found the electric field pro- 
files of the structured mixed phase in cases of spherical and one-dimensional 
geometries. We have demonstrated that the charge screening effect greatly 
modifies the description of the mixed phase, changing its parameters and 
affecting the possibility of its existence. With screening effects included it 
has been proven that the mixed phase structures may exist even if their sizes 
are significantly larger than the Debye screening lengths. The characteristics 
of the structures essentially depend on the value of the ratio of the Debye 
screening lengths in the interior and exterior of the structure and on the 
value of the surface tension on its boundary. We showed the important role 
played by the correlation energy terms omitted in previous works. For real- 
istic values of the parameters these correlation contributions are larger than 
the electromagnetic contribution to the energy discussed in previous works. 

The result, which structure is energetically favorable, essentially depends 
on the values of parameters «o and (3\. This our conclusion essentially differs 
from that was previously done for the Coulomb case, when one ignored the 
screening effects for any j3i and assumed the step-function charge- density 
profile. We recover that limit case result from our general expressions as 
well, but only for tiny values of f3\ and a^ — l. 

We support the conclusion of [16, 17] that the mixed phase may arise 
only if the surface tension is smaller than a critical value. We also evaluated 
the surface tension between the hadron and the quark phases. However this 
evaluation remains rather rough. An uncertainty in the value of the surface 
tension does not allow to conclude whether the mixed H-Q phase exists or 
not. Both possibilities, the mixed phase and the Maxwell construction, can 
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be studied within different models. An advantage of our study is that the 
characteristics of the mixed phase are expressed in terms of only few param- 
eters. Thus knowing these parameters for a given model one may apply our 
results for the discussion of the mixed phase, the Maxwell construction, the 
geometry of the structures, the energy contributions of the structures and 
the electric field profiles. 

In absence of a mixed phase our charged distributions describe the bound- 
ary layer between two separated phases existing within the double-tangent 
(Maxwell) construction. The previously obtained Coulomb limits (when one 
omitted the screening effects) are naturally recovered from our general ex- 
pressions but only in the limiting case of tiny size droplets. This limiting 
case, however, is never realized for realistic values of the parameters. 

In this paper we considered only droplets and slabs and their interplay for 
/ < 0.5 at zero temperature. More detailed analysis of a possible interplay 
between different structures within the whole interval < / < 1 and at finite 
temperatures will be performed in a forthcoming paper. 
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APPENDICES 



A Energy density of quark and hadron sub- 
systems 

In the main text, as a concrete example, we consider the H-Q phase transi- 
tion. We assume that the quark matter (phase I) consists of u, d, s quarks 
and electrons. To easier compare our results for the description of the H-Q 
system with those obtained previously, we choose a simple model of equation 
of state used in [17, 30] with slight modifications. We could use any other 
model. We need it here only to estimate typical values of the parameters. 

The kinetic plus strong interaction energy density of the quark matter 
expressed in terms of particle densities is given by the bag model 



c kin+str 



3tt 2 / 3 



1 + 



2atc 
3tt 



M3 , 4/3 



pT + pT + pT + 



7Y 



4/3 



where B is the bag constant, a c is the QCD coupling constant, m s is the mass 
of the strange quark. Last term is the kinetic energy of electrons. Expansion 
is presented up to the terms linear in a c and m 2 s , where the strange quark 
mass m s (~ 120 -j- 150 MeV) is supposed to be rather small compared to the 
quark chemical potentials, and much smaller values of u and d quark bare 
masses are omitted. 

The hadron matter (phase II) in our model consists of protons, neutrons 
and electrons and the kinetic plus strong interaction energy density is given 
by 

C 37r 2 \4/3 

4in+str " e£nk] + + ^pottPn, P P ] + ™ , (94) 

where e^pi], i — n,p, are standard relativistic kinetic energies of nucleons, 



-kin 



p 2 + m 2 ni(p) 



d 3 p 



n,p, 



(95) 
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m n ~ rrip = itln is the nucleoli mass, rii(p) are step-functional Fermi occu- 
pations. The potential energy e po t we take here in the form 

r I Q (pn-Pp) 2 v . K (p n + p p ) ( p n + p p V 

£pot[Pn, Pp\ = ^0 1" [Pn + Ppjebind H 7^ 1 

Po J-O \ Po / 

+ a a t(Pn+P P )(^^-l), (96) 

Po is the nuclear density (p — 0.16 fm~ 3 ) and values of constants ebmd, K , 
C sat are taken to satisfy the nuclear saturation properties. The value e^nd = 
ebind + ^7V-2po 1 Cp in (PB = Po,P p /Pb = 1/2) provides appropriate binding en- 
ergy §b ind = —15.6 MeV at the nuclear saturation density p . The coefficient 
K is defined by the relation K = K - . 3 f° ■ - 18C sai ~ 285 MeV, where 

_ V m N+PF0 

the value K ~ 240 MeV pro vides an a ppropriate compressibility at satura- 
tion. S = So — 7t 2 po/(4pfo\/pfo + m N) 1S the symmetry energy term. The 
parameter S Q is chosen to be S — 30 MeV, that corresponds to S ~ 18 MeV, 
Pfo = Pf p (Pb = Po,P P /pB = 1/2) = p Fn (pB = Po, P P / Pb = 1/2) is the baryon 
Fermi momentum for the symmetric nuclear matter at the saturation density. 
The value 

C sat = 1p^tf\p B = po, p p /p B = 1/2) - e F0 (97) 



is introduced to reproduce the saturation property, epo = y pfo + m< N- Higher 
order terms ~ (pb/Po — I) 3 i n expansion of the potential energy are for 
simplicity dropped in (96), since their explicit form is not too important for 
our study here. Contribution of mean fields of the heavyjnesons a, uj and p 
is hidden in the values of phenomeno logical parameters K and So- 

Chemical potentials of neutrons and protons are determined with the help 
of (93) - (96): 

i— ; - 2S (p n -pp) K (p n + pp \ 2 

Pn = Pb = yp Fn + mff-\ h e bind + — 1 

Po 6 V po J 

+ K0 ( Pn + Pp _A + 2C sat (Pn+Pp) _ Cgau (9g) 

9 V Po / Po 

Pp = Pn- \fpln + m N + JpFp + m N ~ — ~ ^"(O* 

Po 

p e =(37T 2 Pe ) 1/3 + ^ II (f). 
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As in the phase I, \x e = /x e ,Gibbs for V° = and \i e = for V c 



"/^e.Gibbs- 



Varying these equations of motion, after straightforward calculations we 



get 



C = 



with 



CpnCpp — Ki(C 0p + C 0n ) 
Con — K\ 

45 



0j 



C 0j + 



Po 



(99) 



(100) 



where C 0j - = PFjyPFj + my it is the density of states for the given baryon 
species, j = {n, p}, and 



2S K + 18C sat K oPB 



Po 9p 3pg \ 

Dropping K\ we may use a simplified expression 

Co — Cq p , 



PBy 



(101) 



(102) 



relating to Sp n (f) ~ 0, that is not too bad approximation for p B ~ p and for 
Pb Po, PFn ^> Pf p - We use this simplified expression for rough numerical 
estimates. 



B Surface energy of H-Q interface 

The surface layer separating quark and hadron sub-systems is very difficult 
to describe microscopically mainly due to absence of the microscopic treat- 
ment of the confinement problem. A simplification comes from the fact that 
gradients of all the fields except the electric field are very sharp and this layer 
can be then treated in terms of surface quantities. Thus, in our consideration 
below we will use the fact that in absence of electric field effects the surface 
layer is characterized by a shortest scale ds (^ 1 fm as we estimate). With 
this, we explicitly include electric field and calculate corrections to all the 
meson-nucleon fields due to electric field effects. These corrections, relate 
to larger scales A^ n ~ (5 4- 10) fm corresponding to the Debye screening of 
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the electric field. Due to that A^ ^> d§, the electric field effects do not 
significantly affect the quantities related to the shortest scale ds- 

Let us first discuss the case of the quark-vacuum interface. The quark 
contribution to the energy density of this layer arises due to sharp spatial 
gradients of the quark-gluon fields within the confinement area (5r ~ df ~ 
0.24-0.4 fm) near the bag surface. Since d% is substantially smaller than the 
corresponding length scales A^ n involved in the charge screening problem, 
the quark contribution to the energy of the surface layer can be reduced to 
the surface energy 

E s , q = J dfe s , q = v l <r q \p\ d/R, e s , q = (v l a q d/R)5(fe dD), (103) 

d is the parameter of the space dimension of the structure, d = 3 for spherical 
droplets and d — 1 for slabs, the volume of the quark phase is v l = AttR 3 /3 
for d = 3, and v 1 = 2R for d = 1, and es, q is proportional to the delta- function 
5(f G dD), being zero everywhere outside the surface, / df5(f e dD) = 1. 
The surface tension parameter a q depends on densities of u, d, and s quarks 
in the surface region. 

The surface tension a q [p] is different in dependence of what exterior region 
is considered (vacuum or hadron matter). According to ref. [31], for m s <C 
€p s ~ epn/3 ~ m N /3, where ep s is the Fermi energy of the s-quark, e-p n is the 
Fermi energy of the neutron, the value a q for the quark-vacuum boundary 
(a v q ac ) is estimated as 

< c ^ ^4 s m s . (104) 

/,From eq. (104) we easily evaluate a minimal value of the surface tension 
parameter a q of the quark-vacuum layer as <r mm = m 2 N m s / {12tt 2 ) ~ 0.3m^. ~ 
20 MeV • fm~ 2 (for m s ~ 120 MeV). As one can see, dependence of a q 
on the baryon density is rather weak. For realistic values of densities with 
e Fn ~ 1 4- 1.2 GeV and m s ~ 150 MeV we estimate a q ~ (0.5 4- 0.7)m^. 12 

In the hadronic region gradients of mean fields of the heavy mesons a, uj 
and p are also large compared to gradients of the electric field. As known, 

12 These values are substantially less than the value a q ~ 300 MeV • fm~ 2 used as a 
dimensional estimate in [14] and larger than an estimate a q ~ 10 MeV • fm -2 used in [29] 
in description of the mixed phase structures neglecting the screening effects. 
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thickness of nucleus diffuseness layer is d$ ~ 0.6 fm. This follows both from 
experimental data and from the calculation of the nuclear structure within 
mean field models, cf. [32]. Therefore, contributions arising due to spatial 
gradients of all the hadron mean fields (not related to changes of the electric 
potential, the latter we treat explicitly) can be included in the hadronic part 
of the surface energy density eg,^. This contribution to the surface energy can 
be fixed also with the help of the corresponding surface tension parameter 

E s , h = J dfe s , h = v l a h [p] d/R, e s , h = (v l a h d/R)5(f e 3D). (105) 

The value a^ ac is estimated as ~ 1 MeV • fm~ 2 for symmetric atomic nuclei 
in vacuum. However sharply increases with the density. For the interface 
between two hadronic media a h depends on particle densities in both the 
phases. According to ref. [27], the surface tension between two hadronic 
sub-systems is given by 

ok = a s , h - I e k (re dD-0)- e h (r E 3D + 0)\. (106) 

Here €h\pi\ is the volume part of the hadronic energy density. This expression 
takes into account that the surface tension should disappear (o^ = 0) for the 
case of equal energy densities of the phases. Ref. [27] presented an expression 
for the coefficient a S: h — 0.3p^ 1 (r e 3D — 0)- | psif G 3D — 0) — ps(r e 
3D + 0) | 2 / 3 . As we see, the dependence of as,h on the density is rather weak. 
Therefore, simplifying we will further use constant value as,h- Numerical 
estimation yields as t h — 0A/m n . 

Typical value of the full surface energy per droplet volume v 1 separating 
the hadron and quark phases can be estimated as a half-sum of partial quark 
and hadron contributions, 



E s /v l = j dfe s /v l = dre^/v 1 + J dre^/v 1 



= a[ Pl }d/R, (107) 



where a ~ |(°"g ac + °"ft ac )) °"g ac an d o - ™ are the corresponding values for the 
quark-vacuum and hadron-vacuum boundaries. With above rough estimates 
we would have a ~ (50-=- 150) MeV/fm 2 in dependence on the baryon density. 
This parameter may in general depend on densities of all particle species in 
the surface layer. 
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For a better estimation of surface effects one needs to use chemical equi- 
librium conditions between quark and hadron phases explicitly solving equa- 
tions of motion for the fields in a narrow layer near the bag surface. Avoiding 
this complicated and tedious procedure, we take into account an observation 
of ref. [27] that surface tension should vanish, if energy densities of the phases 
being equal. Thus, analogously to (106) for the phase boundary between two 
quark phases we may assume the validity of the relation 

<y q = a s ,,- | e q (fe dD-0)- e q (f G dD + 0)\, (108) 

where as q is a positive coefficient, which we may estimate with the help of 
(104). 

Unifying results (106), (108) we suggest an interpolation equation 

a = a s -\ e q (fe dD - 0) - e h (f G dD + 0) |, (109) 

with a common coefficient as — |(as,g + a s,h) — (0.3-i-0.4)/m^, and with the 
property that a is zero for e^(r G dD + 0) = e q {f G dD — 0). Here, evaluating 
as we neglected its slight density dependence. 
Thus, we use that 

E s /v I = ad/R, (110) 

with v 1 = AttR 3 /3 for the case of the spherical geometry and v 1 = 2R for 
the flat geometry, and we will estimate the value of the surface tension a 
according to interpolation equation (109). 

Due to dependence of surface tension only on the difference of the bulk 
quantities e q {f G dD — 0) — e/j(r G dD + 0), there appears no contribution 
from the surface region to the chemical potentials (// s = 0). With above 
expressions pressure acquires an additional contribution given by eq. (21) 
but no contribution appears from the density dependence of a. 

Contrary, there is a contribution to the surface charge density 

pf = -de s /dV 

= -asv'd | pl h (dD - 0) - pl(dD + 0) | 5(fe dD)/R. (Ill) 

To get this equation we used eqs (109) and (110). One can see that for 
R 3> l/m n the surface charge, Z s = J dfp 1 ^ 1 , is substantially smaller than 
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all partial contributions to the volume charges given by other terms in r.h.s. 
(17). Integrating (111) yields Z s ~ Z dT (Rm n )- 1 with Z dr =\ Z 1 - Z n |. 
More generally, the typical scale for the change of the electric field is related 
to the Debye screening. Large values of A^ 11 are due to smallness of the 
fine structure coupling constant e 2 entering r.h.s. of the Poisson equation, 
A5 1 ~ em w . Thus, electric field effects may only minor affect quantities, like 
the surface charge, a, etc., relating to a much more narrow surface layer. 
Due to smallness of Zs and to avoid extra uncertainties, we drop the surface 
charge effects, although generalization is straightforward. 

Not entering in more details we permit the variation of the value a in a 
wide region covering the whole interval of parameters used in the literature, 
cf. [17, 9]. 

Concluding, we would like to do a remark. As follows from analysis of 
many works and as was argued above, in the absence of electric field effects 
typical estimate of d$ is < 1 fm. The same estimate is also used in recent 
work [14] studying the first order phase transition from the nuclear matter 
to the color superconductor within the Maxwell construction. On the other 
hand, figures of refs [10, 15] which numerically studied the possibility of 
the mixed phase at the first order phase transition to the kaon condensate 
phase demonstrate a several fermi scale for all the fields: nucleon, sigma, 
omega, rho and kaon fields. There can be two reasons for that. First, these 
works discussing inhomogeneous charged structures within mixed phase do 
not introduce the Poisson equation for the electric field, as well as electric 
field itself, and therefore the region of applicability of their results seems to 
be not quite clear. Second, the kaon field and the electric field in the region 
I relate to the same scale ~ A^, as follows from our analysis of this problem, 
cf. analogous problem of supercharged pion-condensate nuclei discussed in 
ref. [20]. Therefore, surface effects related to the scale ofe, being partially 
smeared by effects related to the larger scales A^ 1 , might be not seen in 
these numerical evaluations. Contrary, a partial smearing of the fields is 
incorporated explicitly in our analysis. Indeed, corrections to all the fields 
due to electric field effects, which we treat explicitly, relate to the scales A^, 11 , 
as we argue, see eqs. (39). In case of the H-Q phase transition that we 
discuss in this paper this smearing can be only partial, since the scale d$ is 
fixed by the quark confinement and can't be substantially changed. 
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C Formulation in terms of Gibbs potential 



In the main text we performed the gauge invariant treatment of the problem 
in terms of the thermodynamic potential and the energy, expressed in terms 
of the pi variables. Sometimes it is more convenient to use another potential 
and variables, e.g. the Gibbs potential G instead of Q. The ground state is 
then determined by the minimization of the Gibbs potential 



where P = const, E is the energy concentrated in the given volume v, [25]. 
With the conservation of the baryon number Nb, appropriate variables of 
the functional G are P, p B , V a , a = I, II, and concentrations of particles. 
Thus, we start with the Gibbs potential functional 



which is expressed in terms of independent variables x u , Xd, V a , p%,x p , where 
x u = Pu/ Pq, Xd = Pd/ p q are concentrations of u, d quarks, baryon density 
Pb in the quark phase p l B = p q /3, p q is the net quark density, x p = p p / p l B 
is the concentration of protons, v 1 is the part of the volume filled by the 
quark phase, v 11 is the part of the volume filled by the hadron phase, v = 
v 1 + v 11 . Densities of strange s quarks and neutrons are fixed by relations 
p s = 3p B — p u — pd, and p n — p l B — p p demonstrating the conservation of the 
baryon number. 

Equations of motion are given by the variation of the Gibbs potential 
functional. Equation 



relates the energy density and its density (volume) derivative to the constant 
pressure P. Other equations of motion are 



G = E + Pv, 



(112) 




(113) 



5G[x u ,x d ,V a ,p a B,x p } 



= 



(114) 



5G[x u ,x d ,V a ,p%,x p ] 
Sx u 



5E[x u ,x d ,V a ,p a B ,x p ] 
Sx u 



= 0, 



(115) 
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5E[x u ,x d ,V a ,p%,x p ] = Q 

= 0, (117) 
= 0. (118) 



5x d 




5E[x u ,x d ,V a 


Pbi x p\ 


5x p 




SE[x u ,x d ,V a 


Pbi x p\ 


SV a 



They do not depend on P since variations are performed at fixed baryon 
density (volume). 

In the presence of the electric field the kinetic and potential contributions 
to the energy density of the quark phase (the energy per the quark fraction 
volume v 1 = f l v ) are given by 

^[x U t x d) V 1 , p l B ] = e[ in+str + e£ inI + €y. (H9) 

The energy density of the quark bag (93) expressed in new variables 
renders 

3tt 2 / 3 / 2a r 



4i„ +str K x d , p\] * B + \- (i + (3^)V3 



X 



4/3 I 4/3 | / -i _ x4/3 _,_ m s(l ~ x u — %d)' 

C « + X d + I 1 X « + 7T 4 / 3 (3p^) 2 / 3 



(120) 



The terms 



= ^+rfCV)W-A(3«.-l)V 1 -^ (121) 

contain an explicit dependence on the electric potential well. To avoid compli- 
cations we selected the way II fixing the gauge (V° = — /i e ,Gibbs) by condition 
V(r -> oo) -> 0. 

The form (121) follows directly from the Lagrangian of the charged par- 
ticles in the electric field. First term in r.h.s. is the kinetic energy density 
of ultrarelativistic electrons, and next three terms are related to the energy 
density of the electric field in the presence of charged sources. The quark con- 
tribution to the electric charge density is p^ 1 = ^p u — \Pd — \Ps = Pb(.3x u — 1). 
The electron density p e is found by a simple counting of all occupied states 
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p e (V) = 2J ^Anp 2 dp/(2TT) 3 ~ -\/ 3 /(37r 2 ), cf. [18]. The last filled level 

(p = Vm = \J (e m — V) 2 — rn%, m e ~ 0.5 MeV is the electron mass) corre- 
sponds to the electron energy e m = m e since particles should not move to 
infinity. For typical values —V ~ m n under consideration a small electron 
mass term can be safely omitted and, therefore, for ultra-relativistic elec- 
trons we get -p e ~ V 3 /3ir 2 . For the case of spatially homogeneous matter 
of a constant density, the electric potential is constant, being determined 
from the local charge-neutrality condition p ch = 0. Then ey = 0, that in 
terms of the electron chemical potential p e ,ioc = —V{f) = const reads as 

e e ,kin = (/i e ,loc) 4 /47T 2 . 

The energy density of the hadron matter (energy per fraction volume 
v 11 = (1 — f l )v) includes several terms 

6 [ X pi Pb\ = e kin+str + e e,kin + e V (122) 

The energy density of the baryon matter contains kinetic and potential con- 
tributions 

e kin+str[ X P> Pb] = 6 kin[ X Pi Pb] + e kin[ X P? Pb] + e pot[^p, Pb], (123) 

now expressed in new variables. The electron kinetic energy density plus the 
electric field energy density, as follows right from the Lagrangian, render 

4i»,ulV"} + 4[x r ,V",p'£} 

= + tionv" - - %£r- (124) 

Equations of motion for the electric field (118) recover the Poisson equa- 
tion (17). Eqs (115), (116) for r G D l + D s yield 

(1 + — ) ( PFu - PFs ) ~ ^ - V + = 0, (125) 

(1 + ^) ( PFd -P FS ) - ^ + ^1^^ = 0, (126) 
3tt 2p Fs 3p B 6x d 

see eqs (119), (120), (121). 

Integrating equation of motion (17) times V{r) and using continuity of 
the electric potential and its derivative at the surface (as above, we neglect a 
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surface charge density), the energy of the quark phase and the surface energy 
given by eqs (119), (120), and (121) can be rewritten in terms of the variables 
Pu, Pd, Ps and p e . Thereby, we reproduce eq. (93). The / dr ^^J term is 
expressed in terms of the particle densities recovering eq. (11). 

With the help of eqs (30) we obtain that equations of motion (116), (117), 
or equivalently (125), (126), specify nothing else than chemical equilibrium 
conditions (32). Surface terms vanish, if eq. (109) is assumed to be fulfilled. 
Then, one can drop surface contributions in chemical equilibrium conditions. 

Integrating equation of motion (17) times V(r) and using the continuity 
of V(f) and W(r) the hadronic part of the energy given by (122), (123), 
(124) is presented in terms of new variables p n , p p , p e , coinciding with (94) - 
(96). The chemical equilibrium condition for the reaction n <-> p + e coincides 
with equation of motion (117) which reads 

JpL^n - + - 2Pp) + V* = 0, (127) 

Po 

where we also used eqs (122) - (124). 

The proton density is not constant in the corresponding eq. (17), rather 
it obeys eq. (127) and, thereby, depends on V(f). In linear approximation in 
deviation of V(r) from its expression for the spatially homogeneous system, 
by variation of (127) we obtain 



c + 45 ° 

Ooj H 



(128) 



Cop + Con] 5 Pp {r) = 5V l \r) + C 0n 5p l *(r), C 0j 

L Po 

Equation of motion (114) acquires the following explicit form 

Pn\PFn + N + Pp V Pf P + 171 N ^ C kin 6 kin 

Po 

+ «(^_ 1 ) + ^ + m + ra = p. (129) 

9 Po \Po J Po 12tt 2 8vre 2 

This equation reproduces nothing else than the standard relation pB-§^ — t = 
P. The same equation of motion is obtained if one takes constant chemical 
potentials from (9) and puts them into equation 

H iPi -e = P. (130) 
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This can be explicitly seen with the help of relations (98) and equation of 
motion (17). Variation of (130) yields 



Using that 



we obtain 



5P = PiSfii + HiSpi - Se. (131) 



5e= {de[pi}/dpi)5pi = piSpi, (132) 



5P = Pl Spi, (133) 



demonstrating that 5P = 0, if 5pi = 0. 

Varying eq. (129) and with the help of (17) expressing the variation 
of the term 5V(f)AV(r)/4iTe 2 (being obtained by partial integration) via 
(p p — p e )V(r), we arrive at the relation 

C 0p Pp - C 0n p n 5p p (r) + C 0n p n - K x p B Sp B (r) - p p 5V(r) = 0. (134) 

/,From eqs (128), (134) we recover the first relation (39). 

We could also vary constant values of chemical potentials (98). Then we 
would obtain 

C JPn(r) - K^pBif) = 0, C 0p 5 Pp (r) - K^Spsif) - 5V(r) = 0. (135) 

Subtraction of these equations reproduces (128). Also with the help of the 
replacement Sp n = 5p B — Sp p we recover (99). 

Finally, we have shown that both methods using Q and G in appropriate 
variables reproduce the very same results. 



D Estimation of non-linear correction terms 

To get a feeling how small might be non-linear corrections to solutions cor- 
responding to the linearized Poisson equation (22), we, following [18], solve 
eq. (17) analytically with a positive charge density p ch = p9(R — r) and 
negative electron density p e = — for simplicity fixing the gauge by the 
condition V(r — > oo) — > 0, for r — ■> oo, p = const ~ p , 6{R — r) = 1 for 
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r < R and zero for r > R. We will also restrict ourselves by consideration of 
the case R ^> Ad, A5 2 = 47re 2 p//x ei Gibbs- Then, in variables x = (r — R)/\^, 
V = —fi e x(x), eq. (17) acquires the form 

X =X 3 -d(-x), (136) 

where reducing 3-dimensional Laplacian to the 1-dimensional Laplacian we 
dropped small 0(Xb/R) term. This equation has exact solutions 

, x < (137) 



X (x) = 1-3 



1 + 



sh(a — x/v3) 
V2 



= J^TVy x > °> ( 138 ) 

with constants a and b found from matching of the potentials with the help 
of the boundary conditions: sha = lly/2, b = 

Linearized eq. (136) expressed in terms of ip ~ 1 — x i n the region x < 
renders 

rj," -3^ = 0. (139) 

Solution is x — 1 ~ Cexp(xy / 3). Using boundary conditions one gets C ~ 
0.24. Comparison of this approximate solution and the exact one (137) 
demonstrates coincidence with deviations < 1.5 %, while the non-linear cor- 
rection 3ip 2 to the l.h.s. of (139) is not small compared to the linear term 
3^. 



E Peculiarities of the Coulomb limit at finite 
quark concentration f l 

In the main text we have demonstrated that our general results coincide 
with those previously obtained in the literature for the limiting case of single 
droplet (/ = f 1 — > 0) of a tiny size. We also have argued that the limit of 
single small transverse size slab is rather peculiar and can't be reproduced 
without inclusion of screening effects. Indeed, the unscreened potential of 
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the single slab diverges at large distances from it. It is however not the 
case for the periodic slab structures at finite /. Let us now consider other 
peculiarities of the case when / is finite, for both geometries. 

Assuming the constant (step function) charge densities p\ h ^ p 1 ^ ref. [16] 
derived the following expression for the energy density of the droplet (d = 3) 
and the slab (d — 1): 



(140) 



Mf) = 



2(d-2) 



-i 



:i __ d/ (i-2/<0 ) + / ( d + 2 )-i. 



Using the global charge neutrality condition (for constant values of p^ h 



i,ii\ 



/Pch + (1 - f)pl = 0, 
we then find for the energies per droplet/slab volume: 



~{d=3) 



47re 2 (pi h ) 2 

5(/-i) 



2 2 . 



(141) 



(142) 



= 2 ^ P Y R2 - (143) 

Adding to these expressions the corresponding surface energy density terms 
floPi/C and p Pi/3£, respectively and minimizing the results with respect to 
£ we arrive at the ratio of the droplet to slab energy densities 

~{d=3) />(d=3)s /, _ 3 f 1/3 I 1 f \l/3 f 1/3 

fc C Igmin J _ q/ 9 /r\l/3 V^ 2^ 1 2 { > { I -i aa\ 

r 11 ^-') ( / ) d-/) 2/3 ' ( ' 

whereas it follows that the ratio becomes to be larger than unit for / > 0.32. 
This means that slabs become to be energetically preferable for / > 0.32 
compared to spherical droplets, cf. the boundary value of the concentration, 
when rods transform to slabs, as it follows from Fig. 17 of [29]. Note that 
again one used here mentioned above "the two part" approach of Glendenning 
to the problem, which allows for the minimization of the energy densities at 
fixed /, see discussion in sect. 7. 
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For spherical droplets (d = 3) and for the constant charge densities we 
obtain 



47reVch r 




Kw \ (145) 



and we immediately recover (142). 

On the other hand, in the case £ <C 1, «i£ <C 1 from our general expres- 
sions (40), (44) we get 

v; ~ ivi'tfc'jv, 

V » = SS^' (r- |?) • (M6) 

Comparing (146) and (145) we express the charged densities through the 
values of the constants V^ 1 and V n that gives the ratio 

Expanding t^ 1 and Vq 1 in (57), (58) in £ <C 1, 5i£ <C 1 and using that 
Vq = Vq, V n = — V n ao£/ f 1 / 3 we reproduce the global charge neutrality 
condition (141) in this case. One can also find the same values of the charged 
densities using expressions (36), (38) in the limit r — * 0, e.g. we get 

(A) 2 = /W(i-/) 2 <»g (147) 

Comparison of (142) with its genetic partner expression (75) in the limit 
1, 5i£ <C 1 also allows to recover the same expression (147). Thus, we 
see that our general equations reproduce the Coulomb limit (145) and the 
expression for the energy (142) for d — 3, for tiny-size droplets. 

As we see from our general solutions, the value p\ h is unambiguously 
determined through the screening parameters k 1 , k 11 and the chemical po- 
tentials, which determine the value of the coefficient f3o, cf. (61). Thus we 
have no room to vary p\ h anymore. In the calculations [16] the value p l (or 
p l ch , which is unambiguously expressed through p l ) is considered as a free 
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parameter, being recovered by the minimization of the total energy density 
(including the Coulomb plus surface energy densities), at the value of the 
concentration /. Since the value p l ch depends on the geometry of the struc- 
ture the value / also depends then on the geometry. However, the latter 
dependence is very weak. Thereby, instead of this more complicated pro- 
cedure with appropriate accuracy one may drop a small contribution of the 
Coulomb and surface energies of the spatial structures to the total energy, 
cf. [8]. Then the quantities p^j 1 and / can be considered as approximately 
independent on the structure geometry, whereas they essentially depend on 
the bulk properties of the matter. In the limit of tiny-size structures p^J 1 and 
/ are given by the same expressions as we get from our general equations. 
Let us consider now slabs. For constant densities and d — 1, we have 

]/ u = ATreyi(xTRw), (148) 

and we immediately recover (143). On the other hand, in the limit £ <C 
l,a 2 t; <C 1 from our general equations (78), (79), (80) we obtain 

v; ~ VoV) 2 *, 

^i-FoV^T^w). (149) 

We suppressed here superscript (1) symbolizing slabs. Comparison of (148) 
and (149) yields 

V?/^ 11 ^ "oPch/Pch- 

Using this relation with the help of eqs (81), (82) we can easily reproduce the 
condition (141). With the help of (83), (84) we can also recover eq. (143). 
Comparison of (143) and eyj + eyjj from (89) leads us to the same expression 
(147) for p l ch , as in three dimensional case. Thus the Coulomb limit is totally 
recovered also in the case d — 1 for finite values of / (£ <C 1, a^C C 1). Other 
limit case, ^< 1, a 2 £ ^> 1, of a tiny fraction volume can't be reproduced by 
the bare Coulomb solution (143), as we have mentioned, leading even to a 
different dependence on £. In this case the correct asymptotic of the 5u)[H is 
given by eq. (88) rather than by eq. (143). 

In our general consideration, as we argued in sect. 7, comparing the 
effective energies 5&iot{£min) f° r droplets and &^tot(£mhi) f° r slabs we recover 
the Coulomb limit only for tiny values of fa. We illustrate it by Fig. 8, 
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where we present the values Sy in + eg 1 " 1 = ey (Cmi n ) + (£min)> as gi yen 
by minimization in £ of the expressions (60), (62), (65), (66) and (69) for 

droplets and £y in + eg 1111 = e^iCmL) + ^(CmhJ) as given by minimization in 
£ of the expressions (83), (84) and (87) for slabs. 
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loc 



Figure 1: Schematic view of the electric potential well with two gauge choices, 
way I: V — V(V° = 0), and way II: V = V(V° = -/x e , Gibbs ); V(way II) = 
V^(way I) - // e ,Gibbs- 
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Figure 2: Fraction dependence of the effective energy per droplet volume as 
a function of the scaled radius £ = R/\\), where (3 is the parameter of the 
energy scale and fii is the parameter proportional to the surface tension, see 
(70), (71). Solid lines are given for / = 0.5 and dashed lines for / = 1/100. 
The Coulomb curves "C" are calculated with the help of eqs (142), (147), 
demonstrating the Coulomb limit, and "e.m." curves relate to "the electric 
field energy" plus surface energy (without "correlation" contributions), a Q = 
Aj} / X% is fixed as «o = 1 . 
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Figure 3: The radius £ m i n = -Rmin/Ap of the spherical droplet corresponding 
to the minimum of the effective energy, 5u3 t ot, in dimensional units. The "C" 
curves are presented for a — 1, as in Fig. 2. 
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Figure 4: The same, as in Fig. 2, but for slabs. The curves labeled by "C" 
are calculated with the help of eq. (88) for / = 0.01 (dashed) and (143), 
(147) for / = 0.5 (solid). 
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Figure 5: The transverse slab radius £ m i n corresponding to the minimum of 
the effective energy, Suj[H, in dimensionless units. The "C" curves are the 
same, as in Fig. 4, shown for ao = 1. 




Figure 6: Comparison of the effective energies of droplets (solid curves) and 
slabs (dashed curves). The numbers at the curves are the values of 
Crossing points of the curves are indicated. 
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Figure 7: The same, as in Fig. 6, but for large values of a . 
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Figure 8: Comparison of the effective energies of droplets (solid curves) and 
slabs (dashed curves) in the Coulomb limit (for (3\ = 10~ 3 ). 
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